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ABSTRACT 

We present Keck high-quality rest-frame ultraviolet (UV) through optical spectra of 21 Type la 
supernovae (SNe la) in the redshift range 0.11 < z < 0.37 and a mean redshift of 0.22 that were 
discovered during the Sloan Digital Sky Survey-II (SDSS-II) SN Survey. Using the broad-band pho- 
tometry of the SDSS survey, we are able to reconstruct the SN host-galaxy spectral energy distributions 
(SEDs), allowing for a correction for the host-galaxy contamination in the SN la spectra. Comparison 
of composite spectra constructed from a subsample of 17 high-quality spectra to those created from 
a low-redshift sample with otherwise similar properties shows that the Keck/SDSS SNe la have, on 
average, extremely similar rest-frame optical spectra but show a UV flux excess. This observation is 
confirmed by comparing synthesized broad-band colors of the individual spectra, showing a difference 
in mean colors at the 2.4-4.4<r level for various UV colors. We further see a slight difference in the 
UV spectral shape between SNe with low-mass and high-mass host galaxies. Additionally, we detect a 
relationship between the flux ratio at 2770 and 2900 A and peak luminosity that differs from that ob- 
served at low redshift. We find that changing the UV SED of an SN la within the observed dispersion 
can change the inferred distance moduli by ~0.1 mag. This effect only occurs when the data probe 
the rest-frame UV. We suggest that this discrepancy could be due to differences in the host-galaxy 
population of the two SN samples or to small-sample statistics. 
Subject headings: supernovae — general, cosmology — observations, distance scale 



1 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA. 

2 Department of Astronomy, University of California, Berke- 
ley, CA 94720-3411, USA. 

Clay Fellow. Electronic address rfoley@cfa.harvard.edu . 

4 Kavli Institute for Cosmological Physics, The University of 
Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA. 

5 Department of Astronomy and Astrophysics, The Univer- 
sity of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637. 

6 Department of Mathematics and Applied Mathematics, Uni- 
versity of Cape Town, Rondebosch 7701, South Africa. 

7 South African Astronomical Observatory, P.O. Box 9, Ob- 
servatory 7935, South Africa. 

8 African Institute for Mathematical Sciences, 6-8 Melrose 
Road, Muizenberg, Cape Town, South Africa. 

9 Center for Particle Astrophysics, Fermi National Accelera- 
tor Laboratory, P.O. Box 500, Batavia, IL 60510, USA. 

10 University of Notre Dame, 225 Nieuwland Science, Notre 
Dame, IN 46556-5670, USA. 

1 1 Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway, NJ 08854, USA. 

12 Institute for Cosmic Ray Research, University of Tokyo, 
5-1-5, Kashiwanoha, Kashiwa, Chiba, 277-8582, Japan. 

13 Institute of Cosmology and Gravitation, Mercantile House, 
Hampshire Terrace, University of Portsmouth, Portsmouth POl 
2EG, UK. 

14 Department of Physics and Astronomy, Johns Hopkins Uni- 
versity, Baltimore, MD 21218, USA. 

lj Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA. 

16 Department of Physics and Astronomy, University of Penn- 
sylvania, 209 South 33rd Street, Philadelphia, PA 19104, USA. 

17 Department of Astronomy and Astrophysics, 525 Davey 
Laboratory, Pennsylvania State University, University Park, PA 
16802, USA. 

18 Institute for Gravitation and the Cosmos, The Pennsylva- 
nia State University, University Park, PA 16802, USA. 

19 Oskar Klein Centre, Department of Astronomy, Stockholm 
University, 106 91 Stockholm, Sweden. 



2 



Foley et al. 



1. INTRODUCTION 

Type la supernovae (SNe la), being both extremely 
luminous and having well-c alibrated peak luminosities 
(starting with the work of iPhillipsI I1993T ) . are excel- 
lent and precise dist ance indicators at cosmological 
scales; see, for example , iLeibundgutl (120011) .iLiviol (2001), 
iPerlmutter fc Schmidt! (|2003t ). and iFilippenkd (|2005bT > 
for reviews. Measu rements of several hundred SNe la at 
low redshifts (e.g . . lHamuv et ah ; 1996; R iess et al.lll999b 
Jha et ail 120061 ; iHicken et al.l l2009at IContrcras ct al. 



2010HGane shalinga m"et al.ll20ldlStritzinger et al.ll201ll) 



and high redshifts have shown that the ex pansion rate 
of the Universe is currently accelerating (jRiess et al.l 
Il998t IPerlmutter et al.l Il999f ). and that the equation- 
of-state parameter of the dark energy, w = P/(pc 2 ), 

is consiste nt with the cosmological constant, w = 

-1 (e.g., iRiess et all I2007T: iWood-Vasey et all [20071: 
Hicken et al.l I2009bt IKessler et al.l I2009at IConlev et al.l 
20111 ; iSullivan et al.l 1201 ll) . The underlying assumption 
of this work is that we can accurately measure distances 
to SNe la at high redshifts using relationships between 
distances and observations of SNe la at low redshifts. 

If there are physical differences between low-redshift 
and high-redshift SNe la (in a way that affects our de- 
rived luminosity distances), we would expect to notice 
discrepancies in their observed spectra. Theoretical stud- 
ies have determined that changing progenitor metallicity, 
a likely difference between low-redshift and high-redshift 
SNe la, produces relatively small variations in the optical 
spectra of SNe la b ut larger differences in their ultraviolet 
(UV) spectra (e g ., lH5fhch et allll998l;ILentz et al.H2000l; 
iSauer et al.l 12001 . IFolev et al.l (|2008aft and lEllis et al l 
(2008) have recently shown that high-redshift SNe la 
are similar to their low-redshift counterparts, though 
both studies had major limitations: an inability to prop- 
erly correct for host- galaxy contaminat ion in the high- 
redshift SN spectra ( Foley et al.ll2008a[) and an insuff i- 
ci ent low-redshif t comp arison sample |Ellis et all feOOS'). 

ISullivan et al.l (|2009) compared composite spectra 
from three samples of SNe la with mean redshifts (z) of 
0.02, 0.48, and 1.16. The spectra in the samples had (by 
construction) similar phases, but the SNe had slightly 
differ ent stretches (con sistent with the trend with red- 
shift; |Howi]r^F|a]l[2007) an d colors. The composite spec- 
tra differed such that the strength of features correspond- 
ing to intermediate- mass elements ( IMEs) are weaker 
at higher redshifts. ISullivan et al.l (|2009D argue that 
the changing stretch distribution with redshift results in 
SN la explosions producing more 56 Ni, and thus hav- 
ing les s mass in IMEs, at higher redshift. ISullivan et al.l 
( 2009) also found that their low-redshift composite spec- 
trum had less flux for A < 3500 A, but they dismiss 
the result since their low-redshift sample contains only 
three SNe l a with spectra covering these wavelengths. 
ICooke et al.l ()201lD created a composite spectrum from 
10 nearby low-redshift SNe la and also found that the 
low-redshift SNe had less flux than the high-redshift SNe 
in the UV. Both the low and high-redshift samples had 
similar light-curve shape distributions, making it unlikely 
that the discrepancy is caused by differences in the pop- 
ulation. Nonetheless, they decided not to examine this 
feature until they obtained more data. 

Further indications of a possible change in the 



UV properties of SN e la come from the analysis of 
IKessler et all (|2009al hereafter K09), who found differ- 
ences in the rest-frame £/-band data for low-redshift and 
high-redshift samples. Using only the Sloan Digital Sky 
Survey-II (SDSS-II) sample to extract cosmological pa- 
rameters, exclusion of the rest-frame UV region caused 
a systematic offset in the measurement of w by ^0.3. 
(However, this offset is significantly reduced when com- 
bining with higher-redshift samples.) This difference is 
referred to as the "[/-band anomaly" bv lK09l . The differ- 
ences could be caused by calibrations (the observer-frame 
U band is notoriously difficult to properly calibrate), or 
by real differences in the samples. Since the inclusion of 
the [/-band and SDSS u-band data had a very different 
effect on the low-redshift sample an d the lowest-redshift 
SDSS-II SNe (z < 0.1), respectively. IK09I suggested that 
the low-redshift, non-SDSS-II sample was the cause of 
the effect, which could easily be the result of incorrect 
c alibration. 

IFolev. Filippenko. fc Jhal (|2008bl hereafter FFJ) 
showed that for a sample of six low-redshift SNe la, the 
steepness of their UV spectra near maximum brightness, 
the "UV ratio" [TZ UV = fx (2770 A)// A (2900 A)], 
appears to be highly correlated with peak luminosity. 
Although this relationship could provide new insights 
into the UV properties of SNe la and reduce the scatter 
in their calibrated absolute luminosities, it is based 
on few SNe (but additional low- redshift SNe l a are 
consistent with the relationship; B ufano et al.l 120091 ; 
IFolev et al.ll2012l ). Moreover, the relationship has not 
yet been tested at high redshifts. 

The SDSS-II SN survey has discovered approxi- 
mately 1000 SNe la at z < 0.4 with ~ 500 be- 
ing spectroscopi c ally confirmed dFrieman et al.l 120081: 
Sako et al.l 120081; IKessler et al.l [2009at iSollerman et al I 



20091: ILampeitl et al.ll2010at ISako et all 1201 ID . We de 

scribe this redshift range as "intermediate," but note 
that if one must choose a binary description, it would 
naturally fall in the "high-redshift" category rather than 
the "low-redshift" one. As a complementary program at 
the Keck Observatory to examine potential differences 
in the UV spectral-energy distributions (SEDs) of low- 
redshift a nd h igh-redshift SNe la and to test the valid- 
ity of the IFFJI relationship at high redshift, we observed 
21 SDSS-II SNe la (the "Keck/SDSS" sample) with the 
intention of obtaining high signal-to-noise ratio (S/N) 
rest-frame UV through optical spectra. We have used 
the deep multi-band photometry of the SN host galax- 
ies from SDSS to correct for host-galaxy contamination. 
Utilizing t he large, high quali ty, low-redshift SN la sam- 
ples fr om IFolev et al.| (|2008af) (desc ribed in greater de- 
tail by IFolev et al.l 12008a! and IFFJI ) , we have been able 
to perform a better comparison of low and high-redshift 
SN la spectra. Not e that the l o w-reds hift composite 
spectra presented by IFolev et ail (|2008aD were created 
using the same me thods as those employed in this paper. 
I Hsiao et al.l (|2007[ ) also constructed a templa te spectrum 
simila r to our low-red shift compos it e spec tra (|Folev et al.1 
l2008al) ; however, the iHsiao et al.l (|2007l ) template spec- 
tra contain high-redshift SN la data and hence are not 
ideal for evolutionary studies. 

This paper is structured in the following way. Section[2] 
discusses our sample of SNe la, the Keck spectroscopic 
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observations, and the data reduction. We also describe 
how we correct for host-galaxy contamination and our 
method for creating composite spectra. In Section [3l 
we compare our Kec k/SPSS sample to the low-redshift 
sample presented by iFolev et alJ (|200 8af) . Section [4] ex- 
plores the implications of our findings for SN cosmology. 
In Section O we split the Keck/SDSS sample by vari- 
ous observational properties to determine their effect on 
spectra. We examine the UV ratio in Section[6] Section[7] 
explores potential systematic errors. We summarize our 
conclusions in Section including some discussion for 
future research. 

2. OBSERVATIONS AND DATA REDUCTION 

The SDSS-II survey was a three-year extension of the 
SDSS (for technic al descriptions o f t he SPSS and its 
data products, see lYork et al.|[2000l and IStoughton et al.1 
l2lffll . The SPSS-II SN Survey, one component of the 
SPSS-II program, was carried out over three seasons, 
each of duration three months (Sep. -Nov. 2005-2007). 
Using the SPSS wide-field imaging c amera (IGunn et ahl 
[19981) on the SPSS 2.5 m telescope (IGunn et alJ 120061) . 
it was a time-domain survey in ugriz over 300 square 
degre es centered on the celestial equator (jFukugita et al.1 
1996). The SN survey took advantage of the excellent 
photo metric and astrome t ric calibration achi e ved by the 
SPSS (IHogg et all f200ll: ISmith et al l [20021: IPier et al.1 
l200l llvezic et al.1 l2004t iTud^" et al. r f2006h . Its main 
purpose was to fill in the "redshift desert" (0.1 < z < 0.3) 
in the existing SN la Hubble diagram (which was vis- 
ible in SN la Hubble dia grams at the time of the cre- 



ation of the project; e.g.. iRiess et al.l [2004J ) . with a fo- 



cus on reducing systematic errors in SN la distance esti- 
mates. In total, the SPSS-II SN Survey discovered ^500 
spect roscopically confirm ed SNe la out to a redshift of 
~0.4 fFricman ct al. 2008?) . The first cosmological results 
were presented by IK09I . with additional inve stigatio ns by 
iSollerman et~aTl (|2009t) and lLampeitl et al.1 (|2010a|) . 

SNe for the Keck/ SPSS sam ple were detected as 
part of the SPSS- II SN survey (jBarentine et al.l 120051 : 
iBassett et al.l2006tl2l)07lfl b1lcT) . In addition to many spec- 
troscopic observations of SNe with smaller telescopes 
flZheng et all 120081 : lOstman et all 120111 : iKonishi et all 
1201 If ), a specific program was developed at the Keck Ob- 
servatory to observe the rest-frame UV spectra of z « 0.3 
SN e la. Targ e ts we re selected in a manner described 
bv iSako et al.l (|2008h : objects which appeared to have 
the appropriate colors and light curves of an SN la near 
maximum light at z « 0.3 were chosen for observation 
in our program. Some host galaxies had predetermined 
rcdshifts, which would occasionally influence the deci- 
sion to observe an object. About half of the targets 
were already spectroscopically confirmed as SNe la near 
maximum brightness in the appropriate redshift range by 
other telescopes before observations were made at Keck. 
No additional data (e.g., host-galaxy properties, bright- 
ness, etc.) were used in deciding which objects were ob- 
served. 

2.1. Photometry 

Precise SPSS-II SN photometry was obtained using 
the "Scene Modeling Photometry" (SMP) technique, 
developed for t h is pur pose and described in detail by 
iHoltzman et~aT1 (j2008D . The basic idea of SMP is to 



model the data as a time-varying point source (the 
SN) and sky background along with a time-independent 
galaxy background, all convolved with a time-varying 
point-spread function (PSF). The model is constrained 
by a global fit that uses all of the images taken during 
the survey that cover the SN region. Each image is mod- 
eled as a sum of fluxes from the SN, host galaxy, sky 
background, and calibration stars near the SN. The cali- 
br ation stars are take n from the SPSS catalog produced 
bv llvezic et al.l (|2007h . The fitted parameters are SN po- 
sition, SN flux for each epoch and passband, and the 
host-galaxy intensity distribution in each passband. The 
galaxy model for each passband is a 20 x 20 grid (with 
a grid scale set by the CCP pixel scale, 0//4 x 0.'.'4) in 
sky coordinates, and each of the 400 x 5 = 2000 intensi- 
ties is an independent fit parameter. In addition to the 
fit parameters, the flux model also depends on measured 
quantities that include sky background, magnitudes of 
calibration stars, position-dependent PSF, and a scaling 
between calibrated flux and CCP counts. There is no 
pixel resampling or image convolution, resulting in cor- 
rect statistical uncertainties. 

2.2. Spectroscopy 

Spectra of SPSS-II targets were obtained with the 
Keck-I 10 m telescop e using the low-r esolution imaging 
spectrometer (LRIS; lOke et al. 1995). Observations in 
2005 were obtained with the 600/4000 grism and the 
400/8500 grating with the polarimetry optics in place. 
Observations in 2006 and 2007 were obtained with the 
400/3400 grism without the polarimetry optics, result- 
ing in better response at blue wavelengths. Additionally, 
an atmospheric dispersion compensator (APC) was in 
place for observations in 2007. This resulted in spectra 
covering the entire optical range, approximately 3400- 
9400 A. All spectra were obtained at relatively low air- 
mass, and several standard stars were observed at similar 
airmass throughout the night. To red uce differential sli t 
losses due to atmospheric dispersion (jFilippenkol 1 1 982ft . 
the spectra were obtained with a position angle near par- 
allactic, except for some objects in 2007 which were ob- 
served at low airmass and with the APC. 

Additional spectra of a subset of the ho st galaxies were 
obtai ned with the MagE spectrograph (Marshall ct al. 
2008) on the Magellan Clay 6.5 m telescope in Pecember 
2008, well after the SNe had faded. The slit was placed 
at the position of the SN with a position angle that went 
through the host nucleus. Because of this alignment, we 
were able to extract spectra at both positions. 

Standard CCP processing and spectrum extraction of 
the LRIS data were accomplished with IRAlPI. T he data 
were extracted using the optimal algorithm of iHornd 
(1986). Low-order polynomial fits to calibration-lamp 
spectra were used to establish the wavelength scale. 
Small adjustments derived from night-sky lines in the 
object frames were applied. For the MagE spectra, the 
sky was subt racted f r om th e images using the method 
described by iKelsonl (|2003ft . We employed IRAF and 

20 IRAF: the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 
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our own routines to flux calibrate the data and re- 
move telluric lines using the well -exposed continua of 
the spectrophotometric standards ()Wade fc HorndH988t 
iFolev et al.ll2003l [2009) . We used the spectrophotornet- 
ric standards BD+284211 and Feige 34 (jOkel I1990T ) to 
calibrate the sp ectra from the blue arm and BD+174708 
and HD 19445 (|Oke fc Gunnlfl98l to calibrate the red 
arm. Our relative photometric flux calibration is typi- 
cally accurate to 5% a cross our entire wavelength range 
(jMatheson et al.l [2008t Silverman et al., submitted). 

The Keck/SDSS sample presented in this paper con- 
sists of 21 SNe la. The properties of the SNe and their 
spectra are listed in Table [1] Two individual spectra are 
presented in Figure [1] (see the online version for all in- 
dividual spectra). Histograms of the sample's redshift, 
rest-frame phase t : and /^Fl are presented in Figure [5J 
It is worth noting that the low and high-redshift light 
curves were fit with different vers ions of the multic olor 
light-curve shape method (MLCS jRiess et aLlll996l ): at 
low redshift by t he publicly available version detailed by 
Uha et al~l (|2007j ). and at high redshift by a slightly mod- 
ified internal SDSS-II version. The differences between 
the implementations are minor, however, and the derived 
values are essentially the same. 

Although all SNe in our sample are obviously SN e la 
based on spectral features (e.g., iFilippenkol 119971 ) . we 
have also perf ormed the Superno v a Iden tification (SNID) 
procedure of iBlondin fc Tonryl (|2007f ) for each spec- 
trum, confirming our classifications. Heliocentric red- 
shifts were determined through either template matching 
(with SNID) or from host-galaxy emission or absorption 
lines, with the precision depending on the method (tem- 
plate matching vs. narrow host-galaxy lines) and the res- 
olution of the spectrum from which the redshift is mea- 
sured. Redshifts determined from template matching, 
from narrow lines in the Keck spectra, and narrow lines 
in either the SDSS main sample or Magellan host-galaxy 
spectra correspond to redshifts with 2, 3, and 4 signifi- 
cant digits in Table Q] and uncertainties of ~0.01, 0.001, 
and 0.0001, respectively. 

The spectra are typically of higher quality than most 
of those previously published for high-redshift SNe la. 
The single exception is SN 16789, which has a relatively 
low S/N and has been removed from our final sample for 
analysis (but its host will be discussed later). There is 
one outlier (compared to the rest of the sample) in both 
redshift and light-curve shape: SN 7147. It is the lowest 
redshift SN (z = 0.111), and separated in redshift space 
from the rest of the sample (the next lowest redshift is 
0.180). It also has a distinct light curve; unlike the rest of 
the sample, which clusters near A = —0.2, SN 7147 has 
A = 0.36, corresponding to a relatively fast-declining 
light curve and indicating a lower-luminosity SN. Its 
spectrum also supports a lower luminosity, with a "silicon 
ratio" of lZ(Si) — 0.4, cor responding to lower - luminosity 
"Branch-nor mal SNe la" (|Nueent et al.lfl995h . 

21 A is a unitless parameter that describes light-curve shape 
and is related to luminosity, with smaller values of A corre- 
sponding higher luminosity; M v (t = 0) = -19.504 + 0.736A + 
0.182A 2 + 51og(Ho/65 km s" 1 Mpc -1 ) mag; A = corresponds 
to Ami 5 (B) = 1.07 mag and stretch s B = 0.96 dRiess et al.lll999 : 
Uha et al.ll2007ft . 
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Fig. 1.— Spectra of SNe 16618 (top) and 7147 (bottom), exam- 
ples of SNe with no and some host-galaxy contamination, respec- 
tively. The black curves are the observed spectra. The blue curve 
for SN 7147 is the reconstructed galaxy spectrum of its host. The 
red curve is the residual spectrum after galaxy-contamination cor- 
rection. The galaxy-subtraction method produces clean SN spec- 
tra. See the electronic version for figures of all spectra of SNe la 
used in this study. 

2.3. Galaxy Subtraction 

As seen in lFoley et al.l ()2008al ). despite careful spectral 
reductions which attempt to remove host-galaxy light, 
many high-redshift SN la spectra have significant host- 
galaxy contamination. Consequently, the intrinsic SN 
SED is difficult to recover from the spectrum. However, 
since the parameter space of galaxy SEDs is well known 
(and well behaved), one can reliably reco nstruct galaxy 
SEDs with broad-band photometry (e.g.. iBlanton et al] 
I2001 . 

Ado pting the approach described by IBlanton et al.l 
(120031) . but upda t ed to include UV wavelengths by 
Blanton & Roweis (2007), and implemented in the IDL 
software package kcorrect . v4_l_4pl. we have used the 
ugriz photometry of the host galaxy at the position 
of the SN, which is a byproduct of the SMP method 
(|HoTt zman et "all 120081 ) , and the redshifts presented in 
Table [T] to reconstruct, using a x 2 minimization method, 

22 http://howdy.physics.nyu.edu/index.php/Kcorrect . 
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TABLE 1 

Supernova Information 
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Redshift 
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(0.02) 


1.1 


52 


2005ik 
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55 
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7 n 


1 i 


zUUojn 


7147 


2U05- 11-05.41S 


0.111 


1.60 


32.3 0.0 


(0.2) 


0.36 


JO. 07) 


0.01 


(0.02) 


12. 


5.6 


zUUojm 


7243 


2005-11-05.319 
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1.25 
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(0.1) 
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JO. 11) 


0.17 


(0.08) 


0.31 


16 


ZUUojl 


1 A V3 

1 4 ( o 
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n oi a 
U.zl4 


1 Q7 
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v u.uy j 
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(0.03) 




( 


2005jn 
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0.33 


1.51 
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-0.16 


J0.08) 


0.01 


(0.01) 






2005jp 
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2005-11-05.544 
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1.64 


32.1 -1.3 


(0.3) 


-0.11 


J0.13) 


0.43 


(0.08) 


13. 


3.6 


2006pr 


16567 


2006-11-23.470 


0.3656 


1.55 


40.4 2.0 


(0.8) 


-0.32 


J0.10) 


0.01 


(0.02) 


12. 


17 


2006pq 


16618 


2006-11-23.447 


0.193 


1.48 


0.0 -4.4 


(0.5) 


-0.36 


J0.13) 


0.08 


(0.07) 


4.0 


51 


2007my 


19027 


2007-10-16.244 
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1.10 


0.0 3.1 


(0.8) 


0.18 


(0.16) 


0.04 


(0.20) 


0.77 


47 


20071u 


19029 


2007-10-14.244 


0.319 


1.13 


25.5 -0.1 


(1.0) 


-0.47 


(0.07) 


0.02 


(0.03) 


2.6 


63 


2007ml 


19101 


2007-10-16.263 


0.19 


1.57 


0.0 -4.0 


(0.3) 


-0.08 


J0.09) 


0.32 


(0.07) 






20071w 


19128 


2007-10-14.267 


0.29 


1.30 


0.0 1.5 


(0.7) 


-0.23 


J0.09) 


0.02 


(0.04) 






2005jo 


7512 


2005-11-05.615 


0.23 


1.89 


-1.0 


(0.4) 


-0.27 


,0.09) 


0.26 


(0.07) 






2006pt 


16644 


2006-11-23.495 


0.2990 


1.49 


-1.9 


(1.3) 


-0.34 


J0.17) 


0.16 


(0.20) 


36. 


8.6 


2006pz 


16789 


2006-11-23.525 


0.3250 


1.54 


5.6 


(1.1) 


-0.36 


(0.17) 


0.15 


(0.15) 


180 


1.1 


2007qu 


20829 


2007-11-12.298 


0.31 


1.14 


-1.2 


(0.4) 


-0.25 


J0.14) 


0.01 


(0.02) 


30. 


2.1 


Note. — SNc 


below the line arc those rejected from the Nominal sample. 


See Scctionl2.4l for the criteria 


for inclusion in 


the Nominal 


sample. 





a Percent galaxy contamination, determined from the synthetic photometry of the spectra in the observer-frame g band. 

Phase relative to the time of S-band maximum light. 
c M v {t = 0) = -19.504 + 0.736A + 0.182 A 2 + 51og(H /65 km s _1 Mpc" 1 ) mag jjha et al.|[2007l) . 

Determined from MLCS fits to our light curves. 
° Determined from kcorrect fits to the host-galaxy broad-band photometry. 

^ Percentage of the total star formation over the last 1 Gyr compared to the total star formation over all time, determined from kcorrect fits to the host-galaxy 
broad-band photometry. 



the galaxy SEDs at the position of the SN. The host 
photometry used to reconstruct the galaxy SEDs is 
therefore derived from the combination of many Stripe 
82 images. Additionally, we em ployed the SDSS DR6 
(jAdelman-McCarthv et all [2008) ugriz photometry for 
the host galaxies to determine galactic properties such 
as star-formation rate and mass. Because of the filter 
gap between the g and r bands, reconstructed SEDs of 
galaxies in the redshift range 0.27 < z < 0.33, in which 
the 4000 A break is between the filters, are slightly de- 
generate. This degeneracy results in many SEDs with 
similar g — r colors but vastly different u — <?, r — i, and 
r — z colors; hence, this method can create unphysical 
SEDs for these galaxies. However, for the typical galaxy 
with z < 0.5, the SEDs are recovered to <0.02 mag i n 
all filters (jBlanton et al.ll2003HBlanton fc Roweisl [2007). 

Since there is more variability in the UV SEDs of galax- 
ies than in optical SEDs, the UV region of the SED 
is more difficult to recon struct than the optical region. 
iBlanton fc Roweil (|2007t ) tested the predictive proper- 
ties of kcorrect by determining the SEDs of galaxies 
with both GALEX and SDSS photometry using only the 
SDSS photometry. They then compared the residual of 
the reconstructed GALEX— SDSS colors compared to the 
observed values. For these galaxies the mean of the resid- 
uals was very close to zero. Since we only wish to recon- 
struct the SED to ~3000 A in the observer frame, we are 
confident that the ugriz data are sufficient for recon- 
structing SEDs over the wavelength region of interest. 
Additionally, since our SNe are found in Stripe 82, our 
host photometry is typically much deeper than that of 
most SDSS galaxies. 

Fortuitously, the hosts of SNe 16644 and 16789 have 



been observed spectroscopically by the SDSS survey 
(jAdelman-McCarthv et al.ll2008| ) , allowing a check of the 
accuracy of this method for reconstructing the galaxy 
spectra. We also observed five host galaxies with the 
Magellan Clay telescope after the SN had faded, where 
we were able to extract a spectrum at the SN position 
and another spectrum at the position of the host nucleus. 
Figure [3] compares the observed and reconstructed spec- 
tra of the hosts of SNe 7147, 7847, 16567, 16644, and 
16789. For the spectra corresponding to the nucleus (in- 
cluding the SDSS spectra), we present the reconstructed 
spectra of the entire galaxy, and not at the position of the 
SN. For the spectra corresponding to the SN position, we 
present the reconstructed spectra at that position. 

The data agree very well over most of the wavelengths 
covered by the observed spectra, suggesting that this pro- 
cedure is a robust way to determine host-galaxy SEDs. 
Since SNe 16644 and 16789 are at z = 0.2988 and 0.3250, 
respectively, we can use these observations to under- 
stand the potential systematic biases for objects with 
0.27 < z < 0.33. For most spectra of both the host 
nucleus and SN position, the observed spectra and re- 
constructed spectra match very well. However, there are 
slight differences for SN 16644. The reconstructed spec- 
trum has more near-UV flux at both the nucleus and 
SN position. Our sample contains eight SNe la with 
0.27 < z < 0.33. However, only five of these SNe are used 
in our final analysis. We have also performed our anal- 
ysis both with and without these SNe (see Section [2~4l) . 
and there is no significant difference in the results. 

After generating a reliable host-galaxy SED with the 
kcorrect routine, we perform two separate methods for 
determining the SN SED. The methods are different in 
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Fig. 2. — Histogram of the redshift (top panel), rest- frame phase 
as determined from the light curves (middle panel), and A distri- 
butions (bottom panel) for the Keck/SDSS SNe la presented in 
this paper. The dashed and open histograms indicate the SNe 
in the Nominal sample and a ddit ional SNe in the All sample, re- 
spectively; see text (Section 12. 4H for a description of these two 
samples. The dashed lines in the bottom panel mark the regions 
of high-luminosity (A < -0.15), normal (- 0.15 < A < 0.3), and 
low-luminosity (A > 0.3) SNe, as defined bv lJha et all 1 120071 ). 



approach, but yield similar results. Both methods are 
described in detail in Appendix |XJ Briefly, "photom- 
etry matching;" i s similar to the method employed by 
lEllis et all (|2008f ). It requires extracting the SN with- 
out attempting to perform local galaxy subtraction and 
warping the resulting spectrum to the combined host 
and SN photometry for an aperture corresponding to 
the combination of slit and seeing. The second method, 
"color matching," is a novel way to remove galaxy con- 
tamination by matching the galaxy-subtracted spectrum 
to the SN colors at the time the spectrum was obtained. 
For our sample, both methods yielded similar results. 
Both methods have their advantages, but the color- 
matching method, which removes as much galaxy con- 
tamination as possible at the time of extraction, should 



et al. 




3000 4000 5000 6000 7000 
Rest Wavelength (A) 

Fig. 3. — Magellan/MagE spectra (black curves) of the nucleus 
of SDSS J232004.44-000320.1 (the host galaxy of SN 7147; top 
panel), SDSS J232004.44-000320.1 at the position of SN 7147 
(second panel), the nucleus of SDSS J020950.32-000342.0 (the 
host galaxy of SN 7847; third panel), SDSS J020950.32-000342.0 
at the position of SN 7847 (fourth panel), the nucleus of 
SDSS J014145.08-001121.8 (the host galaxy of SN 16567; fifth 
panel), SDSS J014145.08-001121.8 at the position of SN 16567 
(sixth panel), the nucleus of SDSS J022716. 17-002336.5 (the host 
galaxy of SN 16644; seventh panel), SDSS J022716.17-002336.5 
at the position of SN 16644 (eighth panel), and the nucleus of 
SDSS J025508.24+001405.2 (the host galaxy of SN 16789; bottom 
panel) . The blue curves are the reconstructed galaxy spectra deter- 
mined by the kcorrect routine. The red curves are SDSS spectra of 
the galaxies. It is worth noting how well the reconstructed spectra 
fit the observed spectra over the majority of optical wavelengths. 
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have a smaller potential for introducing systematic er- 
rors in our sample. We have performed our analysis us- 
ing both methods, and the results do not change signifi- 
cantly. We will thus focus on the color-matched galaxy- 
subtracted spectra. 

2.4. Composite Spectra 

The SN spectra presented in this paper are generally 
of very high quality, allowing for studies of individual 
SNe. However, much additional information can be de- 
termined from a composite spectrum of all SNe. We are 
able to examine both the average properties of the sam- 
ple as well as small-amplitude spectral featur es which are 
below the noise level in individual spectra ()Folev et al.l 
I2008at lElhTeTaI1[2008t ISullivan et al.l[2009h . 

Using the sam e procedure to pro duce the compos- 
ite spectra as in iFolev et aD (|2008a[ ). we created sev- 
eral composite spectra fr om our sample (de tails of the 
sample are presented by IFolev et al.1 l2008al and Silver- 
man et al.; submitted)). This includes removing residual 
emission lines and normalizing the spectra in the region 
4500 < A < 7500 A. The properties of each composite 
spectrum are listed in Table [2] 



The various composite spectra are composed of sub- 
samples of our main sample, culled by quality cuts. 
The "All" spectrum contains the entire sample. All 
other composite spectra exclude the SNe which fail the 
galaxy-contamination correction. The "Nominal" sam- 
ple includes all SNe except those that fail the galaxy- 
contamination routine. Since SN la spectra evolve with 
time, having a small time window will decrease discrep- 
ancies that are the result of large differences in the epoch 
of spectra. Accordingly, we created a "Maximum" com- 
posite spectrum removing the SNe with |t| > 3 days 
(relative to rest-frame B maximum brightness) from the 
Nominal sample. To reduce sensitivity to host-galaxy 
contamination, we created the "Low Galaxy" composite 
spectrum, removing the SNe with high galaxy contami- 
nation (> 8% in the g band) from the Nominal sample. 
Finally, we created a "Low-A" composite spectrum, re- 
moving the low-redshift, high-A SN 7147, which might 
not be representative of the high-redshift population of 
SNe la and could bias our sample, from the Nominal 
sample. 

We have further considered the effect of having SNe 
with 0.27 < z < 0.33 in our sample. Although these 
SNe pass our galaxy-subtraction routines, they may have 
incorrect galaxy subtraction. We have removed these 
SNe from the Low-Galaxy sample to create the "No Gap" 
sample. 

To determine how consistent each composite spectrum 
is with the Nominal composite spectrum in the UV, we 
have determined the ratio of the UV flux in the wave- 
length region 2700-4000 A for each composite spectrum 
relative to the Nominal composite spectrum and listed 
the value in Table [2] The uncertainties given in the 
table are determined from the la boot-strap sampling 
error spectrum relative to the Nominal composite spec- 
trum. These uncertainty measurements do not include 
the boot-strap sampling errors of the Nominal spectrum, 



but we note that the la boot-strap sampling error spec- 
trum for the Nominal sample results in a flux that is 0.94 
or 1.06 (lower and upper bounds) times the flux of the 
Nominal composite spectrum in this region. This metric 
does not contain all pertinent information in the spec- 
trum. For instance, two spectra can have the same av- 
erage flux over a wavelength range but have significantly 
discrepant spectral features over that wavelength range. 
Therefore, we caution use of this number to determine 
the consistency of two composite spectra, but rather to 
use these values as indicators of differences while exam- 
ining the composite spectra (i.e., Figured]) for a detailed 
comparison. 

Figure [4] shows that to within the uncertainties, the 
Nominal composite spectrum is similar to all other com- 
posite spectra mentioned so far, indicating that the cuts 
to the Nominal sample do not affect our results, and 
that the Nominal sample is not significantly biased. Fi- 
nally, the "Strict" composite spectrum excludes all SNe 
in at least one of the following categories: fail the galaxy- 
contamination routine, have \t\ > 3 days, have high 
galaxy contamination, SN 7147, and have 0.27 < z < 
0.33. The Strict sample consists of just three SNe, and 
although they do not have any obvious deficiencies, the 
small sample size means that a single peculiar object may 
significantly bias the sample. 

Examining Figure [4] in detail, we see that the Low- 
Galaxy, No Gap, and Strict composite spectra differ 
slightly from the Nominal composite spectrum. The 
differences between the Low-Galaxy and Nominal com- 
posite spectra are confined to the regions of the spec- 
trum corresponding to the Ca H&K feature, and are 
likely the result of slightly discrepant line velocities or 
line strengths. The differences in line velocities are al- 
most certainl y the result of havi ng slightly different mean 
phases (e.g.- lFolev et al.l 120111 : see also Section l7.1.1|) . 
The UV flux for the Low-Galaxy composite spectrum 
is 1.01 tpio7' further indicating that the differences are 
confined to a spectral feature and do not change the 
continuum. The No Gap and Strict composite spectra 
have more UV flux than the Nominal composite spec- 
trum (having UV flux of 1.20 and 1.22 toll timcs 
that of the Nominal composite spectrum, respectively). 
The No Gap spectrum has slightly more UV flux than 
the Nominal spectrum, but the spectra are consistent 
within the boot-strap sampling errors. The Strict spec- 
trum has a larger deviation from the Nominal spectrum, 
but the Strict spectrum consists of only three SN spectra 
and may be biased. As shown in Section [3l the Nominal 
spectrum has more UV flux than the low-redshift com- 
posite spectrum; the latter has a UV flux 0.84 tom times 
that of the Nominal composite spectrum. The No Gap 
and Strict spectra have larger deviations with the low- 
redshift composite spectrum, indicating that quality cuts 
cannot account for the difference. 

We are confident that the Nominal spectrum is an ex- 
cellent representative of this sample of SNe la. No signifi- 
cant differences are observed when restricting our sample 
to just the SNe near maximum brightness or having low 
galaxy contamination. There is a minor, but insignif- 
icant, difference in the UV continuum of the Nominal 
and No Gap composite spectra. The No Gap sample 
contains only six SNe, and is therefore more suscepti- 
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TABLE 2 

Composite Spectra Information 



Composite No. of 



Supernovae Included in Subsample? 
High- 



Mean 



Mean Mean 



UV Flux 
Relative to 



Name Spectra Rejected 


Contamination 


0.27 < z < 0.33 


|t| > 3 days 


SN 7147 


Phase (days) 


Rcdshift 


A 


Nominal 


All 


21 


Y 


Y 


Y 


Y 


Y 


-0.4 


0.23 


-0.15 


qq +0-07 
u - a3 -0.06 


Nominal 


17 


N 


Y 


Y 


Y 


Y 


-0.3 


0.22 


-0.14 




Maximum 


10 


N 


Y 


Y 


N 


Y 


0.4 


0.22 


-0.12 


1 01 +0 - 09 

J ' ul -0.09 


Low Galaxy 


11 


N 


N 


Y 


Y 


N 


0.0 


0.23 


-0.16 


1 01 + 007 


Low A 


16 


N 


Y 


Y 


Y 


N 


-0.4 


0.23 


-0.18 


1 01 + ' 07 
i - ui -0.06 


No Gap 


6 


N 


N 


N 


Y 


N 


-1.1 


0.19 


-0.17 


1 20 +0 - 05 
liZU -0.27 


Strict 


3 


N 


N 


N 


N 


N 


0.4 


0.20 


-0.16 


-, 99 +0.11 
±,zz -0.14 


Nominal Subsamples 


Active 


7 


N 


Y 


Y 


Y 


N 


-1.2 


0.24 


-0.22 


1 02 +U - UJ 
1 - uz -0.11 


Passive 


6 


N 


Y 


Y 


Y 


Y 


1.0 


0.19 


-0.03 


1 oo +°- 14 

luu -0.15 


Low Mass 


10 


N 


Y 


Y 


Y 


N 


-1.4 


0.23 


-0.18 


n 99 +0 07 
U - 9M -0.07 


High Mass 


5 


N 


Y 


Y 


Y 


Y 


0.9 


0.18 


-0.02 


1 00 + 016 
luu -0.16 


Lower z 


9 


N 


Y 


N 


Y 


Y 


-1.3 


0.19 


-0.12 


1 02 +°- 09 
1 - uz -0.09 


Higher z 


8 


N 


Y 


Y 


Y 


N 


2.2 


0.31 


-0.19 


95 +0- 07 


Negative A 


8 


N 


Y 


Y 


Y 


N 


-1.2 


0.24 


-0.28 


-, n7 +0.08 
-0.10 


Zero A 


9 


N 


Y 


Y 


Y 


Y 


0.4 


0.21 


-0.02 


n 94 +0-10 
U ' M4 -0.08 



ble to biases. Similarly, the Strict composite spectrum, 
which also deviates slightly from the Nominal compos- 
ite spectrum at bluer wavelengths, is derived from only 
three SNe. Given the composite spectra generated from 
various subsamples and quality-cut tests we have per- 
formed, we see no significant biases or systematic errors 
associated with the Nominal sample and its composite 
spectrum. Thus, we will use the Nominal spectrum for 
our analysis. 

3. COMPARISON OF THE KECK/SDSS WITH OTHER 
SUPERNOVA SAMPLES 

3.1. Composite Spectra 

In Figure [5j we present the Nominal composite spec- 
trum and compare it with a low-redshift maximum-light 
composite spectrum. This low-redshift composite spec- 
trum was const ructed from the sa mple of low-redshift 
spectra used bv I Foley et al.l (J2008af). with the ad ditional 
UV spectra of SN 2005cf feufano et al.l 120091) . This 
sample represents one of the largest and cleanest sets 
of low-redshift SN la spectra. Most spectra were ob- 
tained with t he Lick 3 m telesco pe and Kast double 
spectrograph (jMiller fc Stone 1993), which provides well- 
cali brated spectra from the atmospheric limit to 1 /im 
(see lFolev et al.ll2008al and Silverman et al., submitted, 
for additional details). The low-redshift sample used to 
construct the composite spectrum is of very high quality. 
The spectra are of well-monitored (both photometrically 
and spectroscopically) SNe la. Each spectrum is selected 
to have very low galaxy contamination (typically < 3%) 
and very high S/N (typically > 50). The phase and A 
range for the low-redshift sample were chosen to match 
those of the Keck/SDSS sample. To achieve this, the low- 
redshift spectra have —4.8 < t < 3.0 days and A < 0.1. 
From the fourth and fifth panels of Figure El it is clear 
that the average values for phase and A are very similar 
between the two samples at all wavelengths. 

For wavelengths longer than ^5000 A, the Keck/SDSS 
and low-redshift composite spectra are remarkably sim- 
ilar. At wavelengths shorter than ~4400 A, however, 



the spectra differ; the Keck/SDSS spectrum is bluer than 
its low-redshift counterpart. The difference between the 
spectra grows larger with decreasing wavelength, be- 
coming more significant by 3000 A. For wavelengths < 
3400 A, the Keck/SDSS composite spectrum has > 20% 
more flux than the low-redshift composite spectrum. De- 
spite the differences in overall spectral shape, the individ- 
ual features are all consistent between the two composite 
spectra. 

As mentioned above, it is important to note that 
the average properties of the spectra (namely redshift, 
phase, and A) contributing to a given wavelength bin can 
change with wavelength. This is particularly the case for 
the spectra covering slightly different rest-frame wave- 
length ranges. For the Keck/SDSS sample, the change 
of average values with wavelength is the result of different 
redshifts for the SNe, while the average values change in 
the low-redshift sample mainly because of the different 
samples for optical and UV spectra. 

Considering how well matched the Keck/SDSS and 
low-redshift composite spectra are in terms of phase and 
A (the two parameters known to affect an SN spectrum) , 
the differences must come from a different source of vari- 
ation. Specifically, the samples may be mismatched in 
some other parameter not known to cause significant 
color differences, or there is some uncharacterized prob- 
lem with our data reduction techniques. Below we exam- 
ine additional parameters and their effects on the spec- 
tra. 

3.2. Supernova Colors 

The high quality of our spectra presents the oppor- 
tunity to derive properties of the individual SNe. It 
was shown in Section l3~Tl that the Keck/SDSS compos- 
ite spectrum is bluer than the low-redshift composite 
spectrum at UV wavelengths. Although comparison of 
the low-redshift composite spectrum to each individual 
Keck/SDSS SN spectrum is possible, the S/N is much 
lower in the individual spectra. However, by binning the 
spectra, we may be able to sufficiently increase the S/N 
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Fig. 4. — The 1<t boot-strap sampling error region for the Nom- 
inal composite spectrum is displayed in grey. The la boot-strap 
sampling error residuals for the (from top to bottom) All, Maxi- 
mum, Low Galaxy, Low A, No Gap, and Strict composite spectra 
compared to the Nominal composite spectrum are in light blue (X 
minus Nominal). The red curve is the residual of the Lick compos- 
ite and Nominal Keck/SDSS composite spectra (Lick minus Nom- 
inal). The grey region and red curves are the same in each panel. 
The differences between some spectra at ~3650 A and ~6150 A 
are the result of slightly different mean phases, causing a shift in 
the velocity of the Ca H&K and Si II A6355 features, respectively 
(see Section l7. 1,111 , 



of the individual spectra to make an interesting measure- 
ment. 

The binning that we choose is to convolve the spectra 
with various broad-band filters. Figure [6] displays the 
Keck/SDSS composite spectrum with the transmission 
curves of the HST WFC3 F275W and Keck/SDSS ugr 
filters. We also introduce the "NUV" (near-UV) filter, 
a top-hat filter that spans the wavelength range 2700- 
3300 A. The g and r filters probe the area of the spectrum 
that is mostly consistent between the low-redshift and 
Keck/SDSS composite spectra; however, the blue edge 
of the g band does show some differences, so the g — r 
color of the Keck/SDSS composite spectrum should be 
bluer than that of the low-redshift composite spectrum. 
Meanwhile, the F275W, NUV, and u filters probe the 
region where the spectra diverge. 

Using the filter transmission functions, we made color 
measurements for all of the Keck/SDSS spectra. We have 
done the same for the low-redshift composite spectrum. 
Unfortunately, the low-redshift sample rarely has concur- 
rent UV and optical spectra for individual SNe; we there- 
fore are unable to measure the UV minus optical colors 
of individual SNe. To produce reasonable uncertainties 
for the Lick composite spectrum, we performed a Monte 
Carlo analysis of the spectrum using the boot-strap sam- 
pling errors to produce several hundred SN spectra. Our 
measured uncertainties are the standard deviation of the 
colors for these spectra. The color measurements are 
presented in Figure [7] 

Some of our SNe have rest-frame spectra that do not 
extend to sufficiently short wavelengths to cover the en- 
tire filter. We have taken two approaches to measure syn- 
thetic colors for these spectra. First, we set all fluxes for 
wavelengths shorter than the minimum observed wave- 
length in the spectrum to zero. For SNe with spectra 
that extend blueward of 2500 A, where the spectra can 
be particularly noisy, we also set the flux to zero. Since 
any flux at these wavelengths will create a bluer color, 
this measurement sets an upper limit for the color. The 
second approach sets the fluxes for wavelengths shorter 
than the minimum observed wavelength in the spectrum 
to the median flux value for the bluest 100 A of the spec- 
trum. Since the flux generally increases with wavelength 
in the UV, this is likely an overestimate of the flux and 
sets a reasonable lower limit for the color. 

As expected from the comparison of the composite 
spectra, the Keck/SDSS SNe have g — r colors that are 
similar to those of the low-redshift composite spectrum, 
with the composite spectra differing by 0.069±0.025 mag. 
Nearly all of the Keck/SDSS SNe have bluer F275W - g 
and NUV — g colors than the low-redshift composite 
spectrum, with the composite spectra differing by 0.618 
and 0.294 mag for those colors, respectively. Similarly, 
most of the Keck/SDSS SNe have bluer u — g colors 
than the low-redshift composite. Using the Student's t- 
test and assuming that the Lick composite spectrum has 
the correct colors, we find that the probabilities of the 
Keck/SDSS and low-redshift samples having the same 
mean colors are 9. 4xl0 -6 , 1.2xl0~ 4 , and 1.5xl0~ 2 (cor- 
responding to 4.4cr, 3.8cr, and 2.4er) for the F275W - g, 
NUV — g, and u — g colors, respectively. Conversely, 
we find that for the g — r color, the two samples are 
completely consistent with being drawn from the same 
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Fig. 5. — (top panel): Composite spectrum created from our Keck/SDSS sample (black curve) compared to the maximum-light low- 
redshift composite spectrum (blue curve). The spectra are scaled to match in the region 4500 < A < 7500 A. The grey and light-blue 
regions are the la boot-strap sampling errors for the Keck/SDSS and low-redshift composite spectra, respectively, (second panel): The 
grey region is the lcr boot-strap sampling region for the Keck/SDSS composite spectrum. The blue curve is the residual of the Keck/SDSS 
and low-redshift composite spectra. The light-blue region is the residual of the Keck/SDSS composite spectrum and the low-redshift la- 
boot-strap sampling region, (third panel): The number of individual spectra contributing to each wavelength bin in the Keck/SDSS (black 
curve) and low-redshift (blue curve) composite spectra, (fourth panel): The average phase relative to maximum brightness as a function 
of wavelength for the Keck/SDSS (black curve) and low-redshift (blue curve) composite spectra, (fifth panel): The average value of A as 
a function of wavelength for the Keck/SDSS (black curve) and low-redshift (blue curve) composite spectra, (bottom panel): The average 
redshift of the Keck/SDSS composite spectrum as a function of wavelength. 
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Fig. 6. — The composite Keck/SDSS spectrum shown in Figure[5] 
Overplotted are (from blue to red) the transmission curves of the 
HST WFC3 F275W, the newly defined NUV, and the SDSS ugr 
filters. The noise in the Keck/SDSS composite spectrum for rest 
wavelengths shorter than 2500 A is substantial and not plotted. 



population. 

The color information from the individual spectra 
shows that the bluer color of the Keck/SDSS composite 
spectrum is not the result of a minority of SNe heavily 
influencing the overall spectral shape. The probabilities 
also indicate that it is unlikely that the low-redshift and 
Keck/SDSS samples have similar UV colors. 

4. IMPLICATIONS FOR SUPERNOVA COSMOLOGY 

In Section |3j we have shown that the low-redshift and 
Keck/SDSS samples appear to have different UV SEDs. 
We will examine possible explanations for these differ- 
ences, including systematic effects, in Section [7] How- 
ever, for now, we examine some implications of such a 
difference for cosmological results. In particular, we ex- 
plore whether such a mismatch in SEDs could produce 
the [/-band anomaly. 

As discussed in Section [TJ IK09I found differences in 
the rest-frame [/-band data for low-redshift and high- 
redshift samples that, depending on the treatment of the 
rest-frame U band, resulted in a systematic shift in the 
measured value of w. This "[/-band anomaly" is larger 
for some datasets than others. 

We first examine single-band UV distance measure- 
ments in Section I4z.ll finding that a UV excess produces 
a distance bias that is in the opposite direction to that of 
the [/-band anomaly. In Section 321 we find that differ- 
ences between the observed and model UV SEDs signif- 
icantly affect the relationship between light-curve shape 
and observed color. This explain s the relationship be- 
tween the UV flux deficit found bv lK09l and the distance 
modulus offset. 

4.1. Single-Band UV Distance Estimates 

We have shown that there is a mismatch in the UV 
SEDs of the low-redshift and Keck/SDSS samples of 
SNe la. In this section, we show that such a mismatch 
should affect the estimated dista nce to a SN la. However, 
as we will show in Section 14.21 this effect is dominated 
by effects related to the correlation between light-curve 
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Fig. 7. — Color-color diagrams for the individual Keck/SDSS 
SNe la. The black, blue, and red points correspond to SNe with no, 
low (< 8%), and high (> 8%) host-galaxy contamination, respec- 
tively. The black cross with error bars in both directions indicates 
the colors measured for the low-redshift composite spectrum. The 
black "X" indicates the Keck/SDSS composite spectrum. The color 
measurements are derived by convolving the individual Keck/SDSS 
SN spectra with filter transmission curves. Since not all of the spec- 
tra cover the entire wavelength range of each filter, there is a range 
of color for some points. The upper limits are derived by setting all 
flux for wavelengths shorter than the shortest observed wavelength 
of a spectrum to zero. The lower limits are derived by setting all 
flux for wavelengths shorter than the shortest observed wavelength 
of a spectrum to the median flux value of the bluest 100 A of the 
spectrum. 



shape and color when using multiple-band data, and is 
likely not the cause of the [/-band anomaly, but could 
be important for some high-redshift SNe la, which have 
limited color information. 

Some light-curve fitters train on a low-redshift sam- 
ple and apply the measured relationships between light- 
curve shape and color from that sample to high-redshift 
SNe; MLCS2k2 (|Jha et al.ll2007D . in particular, uses this 
methodology. O ther light-curve fitters, such as SALT2 
(|Guv et al. 2003), determine correlations between light- 
curve shape, color, and peak luminosity using all SNe in 
its sample, including both low and high-redshift data. 

MLCS2k2 is trained on the low-redshift sample of 
SNe la. If the average SEDs of low and high-redshift 
SNe la in our observed samples are properly described 
by the spectra presented in Section [3J then MLCS2k2 
will underpredict the flux of a high-redshift SN la in the 
rest-frame UV. 

Imagine that we observed an SN la at high redshift in 
only rest-frame U. The MLCS2k2 estimate of M[/ itrue , 
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the real peak absolute magnitude of the SN in U, will 
be overestimated (MLCS2k2 will underestimate the lu- 
minosity of the SN in U). We can define this difference 
as 

Mtf,true = Mc/,MLCS " AM, (1) 

where M^mlcs is the MLCS2k2 estimate of the true 
value of Mjj,tiue, and AM is the magnitude difference 
between the estimate and the true value. With our as- 
sumption and formalism, AM > mag. 

Since the observed magnitude, U, does not rely on any 
assumptions about the SN SED, we find 

Mtrue = U - M^truc = U- M^MLCS + AM 

= Mmlcs + AM, (2) 

where /Ztruo and a*mlcs are the true and MLCS2k2- 
measured distance moduli, respectively. Since AM > 
mag, MLCS2k2 is underestimating the distance to this 
SN. 

K09 found that MLCS2k2 fits including and excluding 
the rest-frame U band produced a difference in average 
distance moduli for the SDSS-II sample of /i no [/ — A* — 
0.1-0.15 mag, where fJ, no u and fx are the distance moduli 
measured with MLCS2k2 excluding and including the 
rest-frame U band, respectively. If we assume that ig- 
noring the rest-frame U band will avoid the complica- 
tions of any potential SED differences in the UV, then 
fj-noii should be an unbiased measurement of the true 
distance modulus, while fj, represents the measurement 
biased by the rest-frame t/-band data. Under this as- 
sumption, (J, n oU = Mtrue and \x = /xmlcs- Then we have 

AM = /Lttrue — A*MLCS 

~ 0.1-0.15 mag. (3) 

Since AM > mag, the effect of the [/-band anomaly 
is in the same direction as that expected by an excess of 
UV flux in high-redshift SNe la relative to low-redshift 
SNe la. 

To determine if the amplitude of the effect is the same, 
we can simply examine the u — r color of the Keck/SDSS 
and low-redshift composite spectra. For this, we use the 
fact that the SEDs of both composite spectra are very 
similar in the r band; therefore, most differences in this 
color will be from differences in the u band. In Sec- 
tion 13.21 we showed that the difference in u — r colors 
of the composite spectra is A(u — r) = 0.15 ± 0.05 mag. 
This difference is also of the same magnitude as seen by 
iKOl 

Another way in which a difference in UV SEDs can 
affect the measured distance is the treatment of dust. 
Since MLCS2k2 uses the SN colors to estimate the host- 
galaxy extinction, a color difference not accounted for by 
the model will result in an incorrect measurement of this 
extinction. (SALT2 does not attempt to directly measure 
extinction, but the measured distance is still affected by 
SN color in a similar manner.) 

MLCS2k2 models all SNe la with an intrinsic color 
(for our purposes, we will use rest-frame u — r as an 
example) that depends on light-curve shape and has an 
intrinsic scatter. If the measured color is redder than 
the expected color for a given light-curve shape, then 
MLCS2k2 models the difference in color as host-galaxy 



reddening. If high-redshift SNe la have bluer u — r colors 
than low-redshift SNe la (and our measurements of the 
composite spectra support this), then MLCS2k2 will not 
properly measure the host-galaxy reddening, estimating 
a value that is less than the true value. Another way to 
interpret this is that the MLCS2k2 zero-reddening color, 

(u - Oo,MLCS, is t0 ° re d- 

We can describe this effect mathematically. First, we 
define the real extinction in the V band as 

A V ,Uuc = ot[(u - r) - (u - r) ,true] , (4) 

where a is a factor that converts the color difference from 
the zero-reddening value to extinction in the V band (and 
is in part related to Rv), (u — r) is the measured color, 
and (u — r)o.truo is the real zero-reddening color. The 
above effect can be characterized as 

-4y,MLCs = a [(u - r) - (u - r) ,MLCs] 

= a[(u - r) - ((« - r) ,truc + A(u - r))] 
= A v ,true ~ a [A(u - r)] , (5) 

where A(u — r) is the difference between the MLCS2k2 
and real zero-reddening color. Since the real zero- 
reddening color is bluer than the MLCS2k2 zero- 
reddening color, A(u — r) > mag and -Ay,MLCS < 
A v, true- 
Under these simple assumptions, we expect MLCS2k2 
to underestimate the extinction to a given SN la, and 
therefore to underestimate its luminosity and distance. 
Since this effect is in the same direction as the luminosity 
effect described above, the two effects should not cancel 
each other out. 

Although our simple analysis shows that a difference in 
the UV SED for low and high-redshift SNe la, similar to 
what is seen between the low-redshift and Keck/SDSS 
samples, can produce a difference in the inferred dis- 
tances (both in amplitude and direction) as the [/-band 
anomaly, there is a si gnific ant observational problem: 
the SDSS-II sample of K09 has, on average, lower U- 
band flux than their adopted low-redshift sample. This 
is opposite to the difference seen in Figure [5] We dis- 
cuss the potential reasons for the difference between the 
Keck/SDSS and overall SDSS-II samples in Section[7XU 
Below, we discuss how the full SDSS-II sample has the 
opposite UV flux difference (relative to the low-redshift 
sample) yet affects the distance modulus in the same di- 
rection as shown for the single-band case. 

4.2. Multi-Band UV Distance Estimates 

Our presentation of light-curve fitters in the previous 
section is extremely simplified. In reality, all light-curve 
data, including bands other than rest-frame U and data 
obtained away from maximum light, affect the final esti- 
mated distance. Although the difference in the SEDs re- 
sults in a distance modulus difference that is of the same 
magnitude of the U -band anomaly, IK09I calculated dis- 
tance moduli using three filters, and naively, one would 
expect the rest-frame U band to contribute about one 
third of any difference. 

Additionally, colors, which are only measured with 
multi-band data, can influence the inferred distance mod- 
ulus through the correlations between observed color, 
light-curve shape, and peak luminosity. In the MLCS 
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framework, an SN that is bluer than the SN model will 
be assigned a A (the MLCS2k2 light-curve shape pa- 
rameter) that is smaller than it should, and the SN is 
assumed to be intrinsically brighter than it is. That is, 



M true > MmLCS- 



(6) 



From this, we see that the distance modulus is affected 
such that 



Mtru 



m 



M truo 
< m - Mmlcs 

= MMLCS- 



Therefore, 



AM = ^truo - Mmlcs < mag, 



(7) 



(8) 



which is opposite of the effects noted above. 

To assess the impact of these various effects on the in- 
ferred distance, we perfo rmed multiple simula tions using 
the SNANA framework (jKessler et alJl2009bl ). We sim- 
ulated a large sample of UBV light curves for SNe la 
at z = 0.01 with Ay = mag to reduce the affects of 
if-corrections and extinction. 

We fit the light curves in a variety of ways with 
MLCS2k2. First, we fit only the [/-band data. As ex- 
pected (Section 14. X[) . we find AM = /it rU c — Mmlcs > 
mag. Next, we fit the full UBV and only the BV 
light curves with Ay = mag for all SNe. This isolates 
any potential correlations between color-luminosity and 
color-extinction relationships. We find 

Hbv - LHJBV = -0.11 mag, (9) 

which is opposite to the effect described in Section 14.11 
Clearly the multi-band data are important. If we fix 
the t/-band if -correction, Kjju, to zero, the results hold 
to within 0.01 mag; therefore, if-corrections are not a 
dominant factor. 

Next, we fit the light curves fixing the light-curve shape 
parameter, A, for each SN to be the same as the simu- 
lated value for that SN. In this case, we find 



Vbv - P-ubv = +0.05 mag, 



(10) 



in rough agreement with the results of Section 14.11 We 
have performed a similar analysis with an SDSS-like sim- 
ulation (0.1 < z < 0.3 and gri filters) and obtained sim- 
ilar results. 

From these findings, we can conclude that including 
miscalibrated [/-band data in the light-curve fits dramat- 
ically affects the relation between light-curve shape and 
observed color, which in turn affects the measured dis- 
tance modulus. The Keck/SDSS sample has less UV flux 
than the low-redshift sample, on average, and therefore 
has AM > mag. 

5. SUBSAMPLES OF THE KECK/SDSS SAMPLE 

In Scction [2~4l we created several subsamples of the en- 
tire set of Keck/SDSS SNe to test if various quality cuts 
affected the composite spectra. We determined through 
this process that we should exclude SNe where the galaxy 
subtraction procedure failed. 

Since our sample of high-redshift SNe la is relatively 
large, we are also able to create subsamples based on 
physical parameters. As we discussed in Section l3Tl the 
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Fig. 8. — (top panel): Composite spectrum created from 
Keck/SDSS SNe with z < 0.25 (black curve) and z > 0.25 (blue 
curve). The grey and light- blue regions are the lc boot-strap sam- 
pling errors for the lower-redshift and higher-redshift composite 
spectra, respectively, (second panel): The grey region is the la 
boot-strap sampling region for the lower-redshift composite spec- 
trum. The blue curve is the residual of the lower-redshift and 
higher-redshift composite spectra. The light-blue region is the 
residual of the lower-redshift composite spectrum and the higher- 
redshift 1<t boot-strap sampling region, (third panel): The average 
phase relative to maximum brightness as a function of wavelength 
for the lower-redshift (black curve) and higher-redshift (blue curve) 
composite spectra, (bottom panel): The average value of A as 
a function of wavelength for the lower-redshift (black curve) and 
higher-redshift (blue curve) composite spectra. The UV flux is 
similar in both spectra. 



Nominal Keck/SDSS and low-redshift composite spectra 
have very similar average phases and light-curve shapes, 
the dominant factors in spectral differences. Below we 
will separate the Nominal sample based on other physical 
properties to investigate possible physical reasons for the 
difference in UV SED between the Keck/SDSS and low- 
redshift samples. Some properties of the subsamples are 
listed in Table H 

5.1. Redshift 

Splitting the Nominal sample by redshift (using the 
mean redshift of the Nominal sample, 0.25, as the divi- 
sion), we are able to create two samples: a lower-redshift 
sample of nine SNe with (z) = 0.19, (t) = —1.3 days, 
and (A) = —0.12, and a higher-redshift sample of eight 
SNe with (z) = 0.31, (t) = 2.2 days, and (A) = -0.19. 
Composite spectra created from these samples are shown 
in Figured] 

The composite spectra from the lower-redshift and 
higher-redshift subsamples are very similar at most wave- 
lengths. The major difference is in the blue wing of 
the Ca H&K feature, which extends farther to shorter 
wavelengths in the lower-redshift composite spectrum. 
As the velocity of the minimum of the Ca H&K fea- 
ture has a range of about —12,000 to —22,000 km s _1 
at t = -5 days and about -11,000 to -18,000 km s" 1 
at t = 5 days, wi th the line wider at earlier times (e.g., 
iFolev et all 120111 : see also Section 17. 1 . 1|) , the difference 
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in the two composite spectra is likely to be an age effect; 
the lower- redshift sample has an average phase 2.8 days 
earlier than the average phase of the higher- redshift sam- 
ple. 

The other difference between the spectra is at ^4500 A, 
a peak between absorption features. Although it may be 
attributed to the different average redshifts, it is more 
likely that this difference is also an age effect. 

The lower-redshift composite spectrum appears to have 
slightly more UV flux than the higher- redshift composite 
spectrum, but the continua are still consistent within the 
boot-strap sampling errors. Since the SDSS-II sample 
should have relatively little Malmquist-associated biases 
at z fa 0.2, the lack of a difference in UV SED between 
the two composite spectra indicates that these sorts of 
biases are not creating the mismatch in SEDs between 
the low-redshift and Keck/SDSS samples. 

5.2. Host-Galaxy Star Formation and Mass 

Recently, there have been some indications that SNe la 
from different host environments have different peak lu- 
minosities after correc t ion for light-curve s hape and color 
Jfficken et all l2009bl: Kelly et all 120101: iSullivan et ail 
120101 : ILampeitl et all l2010bl ). Although these studies 
have noticed a cor relation between Hu bble residuals and 
both morphology : 1 Ticket i et al.ll2009bf) and mass or star- 
formation rate (Kelly et all 12010 1: ISullivan et al.l 120101: 
ILampeitl et al.l l2010bf) . the cause of the correlation is 
unclear, and it is not known whether the correlation is 
physical or related to the light-curve fitting. 

As a byproduct of fitting galaxy SEDs to the host- 
galaxy photometry, we can obtain a measurement of 
the mass and star-formation history of the host galax- 
ies. The kcorrect routine constructs galaxy SEDs from 
several hundred spectral templates. Using the rela- 
tive importance of each template in the galaxy SED, 
the star-formation history of the galaxy is reconstructed 
(jBlanton fc Roweisl l2007f ) . Once we determine the lu- 
minosity of the galaxy (using the photometry, redshifts, 
and a cosmological model), the spectral templates pro- 
vide a measurement of the stellar mass of the galaxy. 
Although kcorrect does not directly output an uncer- 
tainty f or its mass estimate, the typical uncertainty is 
0.1 dex (jBlanton k Roweisl 120071) . 

From the Nominal sample, there are four SNe with 
no host detected in the SDSS images. Consequently, 
we have no information about the star-formation rate 
of these galaxies (although they likely have high specific 
star- formation rates). We can, with reasonable confi- 
dence, say that the undetected host galaxies are of rel- 
atively low mass. As an example, the host of SN 16567 
at z — 0.37 (the highest redshift galaxy in the sample) 
has a host galaxy with M = 11.8 x 10 9 M Q . Undetected 
galaxies at lower redshift (z = 0.19, 0.215, 0.29, and 0.33) 
almost certainly have lower mass. 

With the full Nominal sample and the subsample of 
13 SNe la having host-galaxy detections, we are able 
to respectively split the SNe by total stellar mass and 
B1000, the percentage of the total star formation over 
the previous 1 Gyr compared to the total star formation 
over all time. This information is listed in Table [T] The 
full sample can easily be divided at M = 5 x 10 9 M© 
into two groups. We include SNe 6933 (z = 0.215) and 




3000 4000 5000 6000 7000 
Rest Wavelength (A) 



Fig. 9. — (top panel): Composite spectra from Keck/SDSS SNe 
with host galaxies that have M < 5 X 10 9 M© (black curve) and 
M > 5 X 10 9 Mq (blue curve). The grey and light-blue regions 
are the la boot-strap sampling errors for the low-mass and high- 
mass composite spectra, respectively, (second panel): The grey 
region is the \o boot-strap sampling region for the low-mass com- 
posite spectrum. The blue curve is the residual of the low-mass and 
high-mass composite spectra. The light-blue region is the residual 
of the low-mass composite spectrum and the high-mass \a boot- 
strap sampling region, (fourth panel): The average phase relative 
to maximum brightness as a function of wavelength for the low- 
mass (black curve) and high-mass (blue curve) composite spectra. 
(fifth panel): The average value of A as a function of wavelength 
for the low-mass (black curve) and high-mass (blue curve) compos- 
ite spectra, (bottom panel): The average redshift of the low- mass 
composite spectrum as a function of wavelength. The two spectra 
are relatively similar, but there is perhaps a slight difference in the 
UV slope, with the high-mass composite spectrum having a flatter 
spectrum from 2900-3500 A. 



19101 (z = 0.19), which do not have detected hosts, in 
the low-mass category. Since SNe 7475 (z = 0.33) and 
19128 (z — 0.29) could conceivably have a mass near our 
chosen division, we do not include either SN in either 
subsample, resulting in a final sample of 15 SNe from 
which we examine the SED dependence on host-galaxy 
mass. The subsample of SNe with B1000 measurements 
is divided into two groups, "active" and "passive" having 
B1000 > 46 and B1000 < 17, respectively. 

Not surprisingly, most passive hosts have a large mass 
and most active hosts have a small mass. Because of the 
larger number of SNe in the mass segregated sample, we 
study that one in more detail. We note that the same 
trends appear when using the passive/active subsamples. 
We present the composite spectra created from the low- 
mass and high-mass subsamples in Figure [§1 

The low-mass and high-mass composite spectra have 
average values of (z) = 0.23, (t) = —1.4 days, (A) = 
-0.18, and (z) = 0.18, (t) = 0.9 days, (A) = -0.02, 
respectively. Although the composite spectra from the 
low-mass and high-mass subsamples (which are essen- 
tially the same as the active and passive subsamples, re- 
spectively) are very similar in the optical, there are some 
minor differences in the UV. In particular, the peak at 
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^3500 A is higher in the low-mass composite spectrum, 
while the peak at ^2900 A is lower in the low-mass com- 
posite spectrum. This results in a relatively flat spectral 
slope for the high-mass composite spectrum in the range 
2900-3500 A. The composite spectra do differ by ^2 days 
in phas e and by ~0-15 in A: however, for the low-redshift 
sample, iFolev et ail ()2008af) did not see the "steepness" 
of the near-UV change as a function of mean phase or 
A (on the other hand, IFFJI did detect a strong corre- 
lation between A and the slope of the SN's spectrum 
between 2770 and 2900 A). We therefore believe that it 
is unlikely that the discrepancies in the mean phase or A 
have caused the differences in these composite spectra. 
Unfortunately, the small sample sizes makes these differ- 
ences marginally significant, and we require more SNe to 
determine whether this trend continues. 

Despite the differences in these composite spectra, 
both have a UV excess compared to the low-redshift com- 
posite spectrum. It is therefore unlikely that the UV 
excess is the result of a biased host population. 

5.3. Light-Curve Shape 

Splitting the Nominal sample by A (using the median 
value of the Nominal sample, A = —0.169, as the divi- 
sion), we are able to create two subsamples, one with 
a smaller average A (the "Negative A") with eight SNe 
and average values of (z) = 0.24, (t) = —1.2 days, and 
(A) = —0.28, and one with a larger average A (the "Zero 
A") with nine SNe and average values of (z) = 0.21, 
(t) = 0.4 days, and (A) = —0.02. Composite spectra 
created from these samples are shown in Figure 1101 

The composite spectra from the zero and negative-A 
subsamples exhibit many differences at all wavelengths. 
The Si II A5972 feature is weak in the negative-A com- 
posite spectrum. The Si II A6355 feature is also sig- 
nificantly weaker in the negative-A composite spectrum 
compared to the zero-A composite spectrum. There are 
also many insignificant differences redward of ^4300 A. 
Clearly the optical SEDs of SNe in these subsamples are 
different, but the individual features are similar enough 
that the spectra are consistent within the boot-strap 
sampling errors. 

The spectra continue to deviate at bluer wavelengths. 
In addition to weak Si II A5972 and Si II A6355 fea- 
tures, the negative-A composite spectrum has a rela- 
tively weak Si II A4130 line. The Ca H&K feature is 
also weaker in the negative-A spectrum. The blue edge 
of the Ca H&K feature has higher flux in the negative- 
A spectrum, indicating that the apparent weakness of 
Ca H&K could be caused by more continuum flux at 
the position of the feature; however, scaling the spectra 
to that peak, the Ca H&K feature is still weaker in the 
negative-A spectrum. Therefore, we believe that this 
difference is inherent to the absorption feature. These 
di fferences can be desc ribed by the lZ(Ca, II) parameter 
of iNugent et al.l (|1995l ). which is the ratio of the flux at 
the red edge of the Ca H&K feature to that of the flux 
at the blue edge. The Keck/SDSS composite spectra fol- 
low the trend of more luminous (or SNe with smaller A) 
corresponding to smaller 7?.(Ca II). These effects have 
also bee n seen when comp aring low-redshift composite 
spectra (jFolev et alj|2008aj ). 

The UV flux is clearly higher in the negative-A spec- 




3000 4000 5000 g 6000 7000 
Rest Wavelength (A) 



Fig. 10. — (top panel): Composite spectra from Keck/SDSS SNe 
with A > -0.169 (black curve) and A < -0.169 (blue curve). 
The grey and light-blue regions are the la boot-strap sampling er- 
rors for the negative-A and zero-A composite spectra, respectively. 
(second panel): The grey region is the la boot-strap sampling re- 
gion for the negative-A composite spectrum. The blue curve is 
the residual of the negative-A and zero-A composite spectra. The 
light-blue region is the residual of the negative-A composite spec- 
trum and the zero-A la boot-strap sampling region, (fourth panel): 
The average phase relative to maximum brightness as a function 
of wavelength for the negative-A (black curve) and zero-A (blue 
curve) composite spectra, (bottom panel): The average rcdshift of 
the negative-A composite spectrum as a function of wavelength. 
Both spectra have similar UV flux. 

trum in the range 3000-3600 A. The absorption min- 
ima of the Fe II A3250 feature is at —23,500 and 
—21,300 km s _1 for the negative and zero-A spectra, 
respectively. This is consistent with the rough trend 
that low-A SNe have a highe r velocity for this feature 
at maximum brightness (|FFJf ). The peaks just blueward 
and redw ard of this feature are called Ai and A2, respec- 
tively, by lEllis et al.l ([20081 . The negative-A composite 
spectrum has a smaller wavelength f or the peak of Ag 
than the zero-A composite spectrum. lEllis et all ([2008D 
found that there was significant scatter in the position of 
these features and that the position was strongly corre- 
lated with phase. 

Examining Figure 15 of lEllis et all §008), we see that 
the Keck/SDSS sample shows similar results when sep- 
arating the sample into two bins based on light-curve 
shape. Their "high-stretch" (corresponding to smaller 
A) spectrum has stronger UV flux, a smaller 7^(Ca II), a 
higher velocity for the Fe II A3250 feature, and a shorter 
peak wavelength for A2. All of these characteristics are 
also seen in the Keck/SDSS sample. 

6. INVESTIGATION OF THE UV RATIO 

The set of Keck/SDSS SN la spectra is comparable 
to t he cu rrent low-redshift UV sample in size and qual- 
ity (jFFJT ). The low-redshift UV sample was mostly ob- 
tained with the International Ultraviolet Explorer (IUE) 
telescope, including almost all spectra n ear maximum 
light (SN 2005cf was observed by Swift; iBufano et al.l 
2009). There are 15 low-redshift SN la UV spectra with 
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Fig. 11.— Rest-frame UV spectrum of SN 16618 (z = 0.193; 
t = —4.4 days. Marked are the wavelengths (2770 and 2900 A) 
and corresponding fluxes used for calculating the UV ratio, TZyy 
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Fig. 12. — Relationship between the UV ratio (R.jjv) an d A. The 
black points represent the SNe in our Keck/SDSS sample. The blue 
squares are the low-redshift sample from FFJJ. The Keck/SDSS SNe 
are obviously offset to higher values of TZjjy, the result of the UV 
excess. 



-3 < t < 3 days from 7 SNe la (SNe 1980N, 1981B, 
1986G, 1990N, 1991T, 1992 A, 2005cf). These "classic" 
SNe la all have precise distance measurements from inde- 
pendent Cepheid observa tions or surface brightness fluc- 
tuations. Consequently, iFFJl were able to determine a 
relationship between the ratio of the flux of the SN spec- 
trum at 2770 and 2900 A, TZ UV (the "UV ratio"), and 
peak luminosity. An example of how the UV ratio is 
calculated is shown in Figure [TT1 

The UV ratio is an excellent luminosity indicator for 
these six low-redshift SNe la. Using a linear relationship, 
we see that the UV ratio yields a scatter in the peak ab- 
solute magnitudes of the SNe of 0.21 mag, comparable 
to the scatter found when using the luminosity vs. light- 
curve shape relationship (a « 0.17 for this sample; iFFJf ). 
Furthermore, combining the UV ratio and the luminos- 
ity vs. light-curve shape relationship, IFFJl were able to 
reduce the scatter to 0.09 mag. 

For a subsample of our Keck/SDSS SNe, we possess 
accurately fiuxed-calibrated UV spectra (based on our 
galaxy-subtraction techniques) that have reasonable S/N 
at 2770 A to measure TZjjy. These SNe are all SNe from 
the Nominal sample observed in 2006 and 2007, supple- 
mented by those with z > 0.215 from 2005. In Figure [T2l 
we present the measured UV ratios from our sample as 
a function of A and compare to the low-redshift values. 
The SDSS-II light curves were fit with a slightly differ- 
ent version of MLCS2k2 than that of the low-redshift 
sample (see Section |2~2"1) . Fitting a subset of low-redshift 
SNe with both versions of MLCS2k2, we see an offset 
of 0.01 and a root-mean-square scatter of 0.02 in A be- 
tween the two outputs. The larger uncertainty in A for 
the Keck/SDSS SNe is the result of slightly lower qual- 
ity (both S/N and temporal sampling) light curves. The 
larger uncertainty in TZjjv for the Keck/SDSS SNe is 
due to the combination of slightly lower-S /N spectra and 
larger values of TZjjv (for spectra with a fixed S/N, one 
naturally expects a larger uncertainty for larger values of 
K uv ). _ 

Examining Figure 1121 we immediately notice that for 



a given A, the Keck/SDSS sample has a higher value 
of TZjjv- This is expected in view of the UV excess in 
the SDSS-II sample noted above; the UV excess is more 
prominent at shorter wavelengths. Unfortunately, the 
small range of A for the Keck/SDSS sample does not 
allow us to determine whether the trend seen at low red- 
shift (i.e., TZjjv increases with A) also exists at high red- 
shift. 

As a test of using TZjjv a s a luminosity indicator, we 
wish to determine My for each Keck/SDSS SN la and 
compare that to TZjjy . To find My, we first measure 
the peak apparent r or i magnitude (depending on the 
redshift of the SN) using MLCS template light curves 
fit to the photometry in that single band. Then, using 
the concordance cosmology (Hq — 70 km s -1 Mpc -1 , 
Cl m = 0.3, 51a = 0-7, and w = —1), we can determine 
the distance modulus ( u) of each SN. From our MLCS 
fits and iSchleeel et al.l (| 19981) . we are able to determine 
the host-galaxy and Milky Way reddenings, ^^host and 
Ax.mWj respectively. Finally, from our MLCS fits, we 
are also able to derive a if-correction, Kxv, where XV 
denotes a correction from observed-band X to rest-frame 
V. Then, using the equation 



My = X n 



.4 



V.host 



.4 



v.mw 



K 



XV i 



(11) 



where X denotes either the r or i band, we are able to 
determine the absolute visual magnitude for each SN. 
Note that there are assumptions made regarding the K- 
corrections which may not be appropriate for our sample 
given the differences between the spectra of low-redshift 
and Keck/SDSS SNe seen above. However, the rest- 
frame V band is relatively free of these differences. An- 
other assumption is that the MLCS-measured value for 
Ay applies to our sample given the observed color differ- 
ences between the high and low-redshift samples. How- 
ever, uncertainties in Ay and K r y are typically ~0.05 
and 0.02 mag, respectively, which are generally smaller 
than the errors associated with the determination of r max 
and \i. 
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TABLE 3 
My — IZuv Linear Fits 



Fig. 13. — Relationship between the UV ratio (IZuv) an d My 
at peak brightness. The black points represent the SNe in our 
Keck/SDSS sample. The blue squares arc the low-redshift sample 
from lFFJl There are multiple measurements for some low-redshift 
SNe (from different spectra over — 3 < t < 3 days). Each low- 
redshift SN has a unique value of My . The solid black line is 
the relationship given in Equation 1131 The dashed black line is 
the relationship given in Equation 1141 including the low-redshift 
data with My < —19.0 mag, corresponding to smaller A. The 
Keck/SDSS SNe follow a different trend than the low-redshift SNe 
that cannot be attributed to only the UV excess, which would 
shift the SNe to higher values of TZuv (and potentially change the 
slope slightly). However, it is possible that all SNe with My < 
— 19.0 mag follow a simple relationship. 
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(b) the low and high-redshift samples have truly different 
relationships, and (c) subsamples of SNe form groups in 
the TZuy-My plane. 

7. DISCUSSION 

To summarize our results, we find (Section [3]) that 
the intermediate-redshift Keck/SDSS sample has signifi- 
cantly more flux (> 20%) in the near- UV (A < 3400 A) 
relative to a low-redshift sample. This mismatch has sig- 
nificant implications for measuring distances with SNe la 
(Section [4]). In Section [5j we examine subsamples of the 
Keck/SDSS sample to determine if there were any obvi- 
ous correlations between spectral features and rcdshift, 
host properties, and light-curve shape, respectively. We 
find that the UV SED of the Keck/SDSS sample does not 
depend on redshift, and there are no significant differ- 
ences that correlate with host-galaxy mass or light-curve 
shape. 

We now investigate possible systematic errors which 
could create a UV mismatch. We also discuss this result 
in the context of SN evolution. 



We compare the values of TZjjv an d My for the 
Ke ck/S DSS and low-redshift SNe in Figure US As found 
by IFFJI . the low-redshift SNe follow the relationship 

M v (t = 0) = -19.366 + h.Qm{K uv ~ 0.3) mag, (12) 

where we have adjusted My for the SNe to match H = 
70 km s -1 Mpc -1 . It is obvious that the trend seen in 
the low-redshift sample is not present in the Keck/SDSS 
sample. One might expect a shift of the trend to higher 
TZuv, but a similar slope, if the UV excess is equal for 
all Keck/SDSS SNe. However, fitting all Keck/SDSS SNe 
with spectra that reach rest-frame 2750 A, we find 

M v (t = 0) = -19.702 + 0A35{K UV - 0.3) mag, (13) 

which is much shallower than the low-redshift trend. The 
fit corresponding to Equation [T3] has x 2 /dof = 15.9/11. 

The low-redshift sample contains only six SNe, while 
the Keck/SDSS sample consists of SNe from a small 
range of My (and A). Four of the low-redshift SNe ap- 
pear to be consistent with the relationship derived from 
the Keck/SDSS SNe, while two (those with the highest 
value of A, and corresponding to the five low-redshift 
points with My > —19.1 mag in Figure [TBI are incon- 
sistent. If we examine all SNe (at both low and high 
redshift) with My < —19.1 mag, we determine 

M v (t = 0) = -19.546 - 0.110(Kuv ~ 0.3) mag, (14) 

with x 2 /dof = 25.9/18. The values for the linear fits 
for Equations I12H14I and corresponding uncertainties are 
presented in Table [31 

Both samples must be expanded to determine whether 
(a) these trends are real within their own populations, 



7.1. Possible Systematic Errors 

There are many possible systematic errors associated 
with the spectral shape of the SNe la in our study. They 
can be separated into three categories: sample selection, 
spectral reductions, and post-processing techniques. We 
address each of these in turn. 

7.1.1. Sample Selection 

The above results rely on the assumption that the 
Keck/SDSS and low-redshift samples were randomly se- 
lected from their respective populations. Here we discuss 
a possible source of systematic uncertainty from selection 
bias. If the Keck/SDSS sample presented here is not rep- 
resentative of high-redshift SNe la, then any conclusions 
drawn from this sample are possibly incorrect. Similarly, 
if the sample used to create the low-redshift composite 
spectrum is biased relative to the sample of SNe used in 
a cosmological analysis, the same caveats apply. 

If the Keck/SDSS sample is not a random subset of 
the larger spectroscopic SDSS-II sample, then we would 
expect there to be some additional observational evi- 
dence besides a UV excess. In Figure HU we compare 
the light-curve shape parameters of our subsample to the 
99 SN e la used in the first SDSS-II cosmological analysis 
( K09). Performing a Kolmogoroff-Smirnov (K-S) test, we 
find that the two samples are not significantly different, 
with a p-value of 0.348. Therefore, the Keck/SDSS sam- 
ple appears to be representative (in terms of light-curve 
shape) of the overall SDSS-II spectroscopic sample. 

In Figure [TBI we show the measured value of A for each 
SN as a function of redshift. There is no real indication 
of Malmquist bias. A best-fit line to the A distribution 
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Fig. 14.— Histogram of the A distribution of the SDSS-II SNe la 
from the first SDSS-II cosmological analysis and the subsample 
presented in this paper (dashed histogram). The dashed lines mark 
the regions of overluminous (A < —0.15), normal (—0. 15 < A < 
0.3), a nd underluminous (A > 0.3) SNe, as defined by IJha et al.1 
(2007) and shown in Figure [2] 
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Fig. 15.— The A distribution of Keck/SDSS SNe la presented in 
this paper as a function of redshift. The filled black diamonds are 
the Keck/SDSS SNe in the Nominal sample. The open black dia- 
monds are the Keck/SDSS SNe rejected from the Nominal sample. 
The smaller blue diamonds are the full SDSS-II first-year sample. 
A K-S test shows that the Keck/SDSS and full SDSS-II samples 
are consistent. 



of the Nominal sample with redshift is l.ler from zero 
slope. This, combined with the changing demographics 
of the SN p opulation and the effect on average light- 
curve shape ([Howell et al.l 12007). suggests that there is 
no significant Malmquist bias. 

The full, photometric SDSS-II sample is certainly 
Malmquist biased, but the SDSS-II sample should be 
nearly unbiased to z ~ 0.3 for discoveries (as opposed 
to sp ectro scopic completeness, which is limited to z w 
0.15; IK09D ; the Keck/SDSS sample may be similarly 
Malmquist biased, although we do not strongly detect 
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a Approximate rest-frame wavelength of the center of 
the g band at the average redshift for the redshift bin. 



it. 

Besides a simple Malmquist bias where we are select- 
ing more-luminous SNe, we may be selecting bluer SNe 
than the typical high-redshift SN la included in cosmo- 
logical analyses. This is particularly problematic for 
z > 1 SNe la where searches are usually performed in 
the rest-frame near-UV (but photometrically followed in 
redder bands). For the SDSS-II sample, the search is 
performed in gri, avoiding the near-UV, which should 
reduce such a bias. The difference between g — r color 
of the Keck/SDSS and low-redshift composite spectra is 
A(g —r) = 0.1 mag at z — 0.35, and less for lower red- 
shifts. Therefore, it is highly unlikely that such a bias is 
directly affecting our sample. 

Nevertheless, the Keck/SDSS sample has an excess of 
rest-frame UV flux relative to the low-redshift sample, 
which in turn has an excess of re st-frame UV flux relative 
to the full SDSS-II sample (K09). Examining the g-band 
MLCS2k2 light-curve fits for the lK09l sample, we can as- 
sign a value, which we call puv, of —1, 0, or 1 to light 
curves that are below, consistent with, and above the 
MLCS2k2 g-band model, respectively. Our results are 
shown in Table |4] With increasing redshift, correspond- 
ing to bluer wavelengths, the g-band data of the SDSS-II 
sample go from being above the model to significantly 
below the model. The 1 2 SN e from the Keck/SDSS sam- 
ple that are also in the IK09I sample, on the other hand, 
have (puv) = 0.23 (indicating that the sample tends to 
be above the model), with (z) = 0.250. For the 0.2- 
0.3 redshift bin, the IK09I has (z) = 0.25 (the same 
as the Keck/SDSS sample within the IK09I sample) . but 
(Puv) = —0.31. Although there is no obvious color se- 
lection bias that causes this effect, the differences appear 
to be real. It is therefore important to examine other 
differences between the various samples. 

The low-redshift sample has few UV spectra, and it 
is possible that these SNe are not representative of low- 
redshift SNe la. The UV portion of the maximum-light 
comp osite spectrum consists of 15 spectra from 7 SNe 
(|FFJ| : SNe 1980N, 1981B, 1986G, 1990N, 1991T, 1992A, 
and 2005cf), all observed with IUE or Swift. The low- 
redshift composite spectrum begins to deviate from the 
Keck/SDSS composite spectrum at ^4000 A, which is 
easily observed from the ground (and not included in any 
IUE spectra), is not heavily affected by the atmosphere, 
and is sampled by many spectra. We therefore do not 
think that the small sample of low-redshift UV spectra 
is biasing our results. Although a larger number of low- 
redshift UV spectra of SNe la is desperately needed (and 
is being addressed by our Swift program^ and comple- 

23 PI Filippenko, Programs 04047, 5080130; PI Foley, Program 
6090689. 
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mentary HST programs^, a biased low-redshift set of 
SNe la having available UV spectra does not appear to 
be the cause of the UV difference. 

Alternatively, the low-redshift and Keck/SDSS sam- 
ples might have intrinsic differences that may cause dis- 
crepancies in UV spectra but not show differences in 
light-curve shape. To properly compare the samples, 
we must consider the intrinsic differences of the surveys 
that discover SNe at low redshift compared to that of 
SDSS-II. The predominant method of finding SNe at low 
redshift is through targeted surveys, where a sample of 
galaxies is searched continually. In particular, this is the 
method used by the Lick Obse r vatory Supernova Search 
(LOSS: iFihppenko et al.l[200lt lFilippenkoll2005a)) . from 
which a large fraction of the low-redshift SNe la come. 
SDSS-II, on the other hand, is an untargeted survey, ob- 
serving the same part of the sky continually, but not 
targeting particular galaxies. This means that the sam- 
ple demographics are different. Although most of the 
stars in an untargeted survey are in galaxies similar to 
those searched by targeted surveys, so they are not dra- 
matically different, there may still be a significant differ- 
ence in the host-galaxy population. Specifically, the host 
galaxies of the Keck/SDSS sample have a median mass 
of (2.6-3.6) xlO 9 M (depending on how one treats the 
host- galaxy nondetections), while the average masses for 
E, Sb, and Irr galaxies from LOSS, from which a signif- 
icant number of the low-r edshift sample come , are 16.0, 
6.4, and 2.4 x 10 10 M Q (|Leaman et al.ll2010D . There- 
fore, the average SN la from the Keck/SDSS sample 
is hosted in a galaxy that is probably ^0.1 times the 
mass of the average irregular galaxy and ~0.02 times the 
mass of the average elliptical galaxy in the low-redshift 
sample. Given the mild differences found in the spec- 
tra of SNe hosted in low and high-mass hosts (see Sec- 
tion [OJ and the corr elation between host-galaxy mass 
and Hubble residuals (iKellv et al.l 120101 : iLampeitl et al.l 
120101* ISullivan et all 120101 ). we consider variations im- 
printed on SNe la by their environments as a possible 
cause of the difference. However, the Keck/SDSS sample 
is a subset of the SDSS-II sample, and it is unclear how 
it could have significantly different average host proper- 
ties. Future untargeted, low-redshift SN samples from 
surveys such as the Palomar Transient Factory (PTF), 
La Silla-QUEST, SkyMapper, or Pan-STARRS should 
provide SNe with similar host-galaxy properties as the 
high-redshift SNe . (However, the low-redshift sample of 
ICooke et al.ll20lil primarily came from PTF and showed 
similar UV flux as our low-redshift sample.) 

A second difference between the searches is that those 
at low redshift tend to be done in a single band (of- 
ten unfiltered, which for LOSS is similar to the R band), 
while the SDSS-II search is performed in gri. The SDSS- 
II search requires a detection in two of the filters. At 
z « 0.25, these are typically r and i. It is therefore un- 
likely that the Keck/SDSS sample is biased to select sig- 
nificantly bluer SNe at discovery than for the low-redshift 
sample. In a typical exposure, the SDSS-II SN survey 
achieved a 50% detection limit of g — 22.7 mag, which 
is ^0.5 mag deeper than the typical z = 0.35 SN la at 
peak. As mentioned above, the difference between the 
g — r color for the Keck/SDSS and low-redshift compos- 

24 PI Ellis, Programs 11721 and 12298; PI Foley, Program 12592. 



ite magnitudes is A(g — r) — 0.10 mag at z = 0.35 (and 
A(g — r) < 0.2 mag until z w 0.55), so we should not be 
significantly affected by this potential bias. Additionally, 
the rest-frame g—r colors of the Nominal and low-redshift 
composite spectra is only 0.069 mag. 

Examining FigureO we see that the Keck/SDSS sam- 
ple is predominantly composed of SNe with broad light 
curves (i.e., low A). Although the full low-redshift and 
Keck/SDSS samples have different demographics, it is 
important to point out that (a) this is a natural conse- 
quence of t he evolving demogr aphics of SNe la with red- 
shift (e.g., [Howell et all 120071 ) . and (b) the low-redshift 
sample we have used to generate the comparison low- 
redshift composite spectrum has very similar properties 
to that of the Keck/SDSS sample (as demonstrated by 
the comparison of phase and A in Figure [5]). Conse- 
quently, differences in the spectra as a result of different 
photometric and phase properties should be largely mit- 
igated and not biasing our findings. However, note that 
we draw some conclusions in later sections based on the 
full low-rcdshift photometric sample which may not show 
a similar trend. 

7.1.2. Spectral Reductions 

Another set of possible systematic errors is associated 
with the observations and reduction process. The two 
major sources of systematic error would be the system- 
atic "reddening" of the spectra as the result of the spec- 
troscopic slit not b eing aligned along the parallactic angle 
(Filippcnkol ll982| ) and incorrect flux calibration. 

All of the Keck observations were taken within 5° of 
the parallactic angle to preserve the intrinsic SED of the 
SN spectrum. Furthermore, the observations taken in 
2007 were with the ADC in place, reducing the effec- 
tive atmospheric dispersion to a negligible amount for all 
but the highest airmasses. Moreover, if significant atmo- 
spheric dispersion were present at high airmass, the blue 
flux should decrease more than the red, given that the 
SNe were centered in the slit while using the red-sensitive 
LRIS guide camera — yet this is the opposite of the ob- 
served difference between the low and high-redshift sam- 
ples. (Note that the standard stars were observed in the 
same manner as the SNe.) In addition, the low-redshift 
UV spectra were obtained from space, where atmospheric 
dispersion does not affect the observations. 

The second possibility is that the flux calibration of 
the SN spectra was incorrect, causing an overestimate 
of the bluest fluxes. However, the spectra were obtained 
over three seasons on several nights. Therefore, any flux- 
calibration error must have occurred for these tempo- 
rally separated observations. The Keck/SDSS SN spec- 
tra were also reduced in exactly the same manner as 
that of the optical spectra from the low-redshift sam- 
ple (and similar in methodology for the UV spectra), 
with all Keck/SDSS SN spectra and a significant number 
of low-redshift SN spectra reduced by the same person 
(R.J.F.). If there is an inherent flux-calibration error for 
all spectra, this would have been observed in the low- 
rcdshift sample as well. A second set of spectral reduc- 
tions was performed by another person (Ryan Chornock) 
with similar results. Additionally, the spectra obtained 
in 2005 were observed with polarimetry optics, produc- 
ing two spectra for each object that were reduced with 
completely separate calibrations, yet they always agreed 
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within the uncertainties. A final argument against a flux- 
calibration error is that these errors would be present in 
the observer frame, not in the rest frame where the differ- 
ence is apparent. Since our spectra span a range in red- 
shift, the errors should be present at varying wavelengths 
and not correspond to the same rest-frame wavelengths. 
However, when separating our sample by redshift, we see 
no difference (see Section I5.1|) . Because of our metic- 
ulous flux calibration with multiple decades of testing, 
flux-calibration errors seem unlikely. 

7.1.3. Post-Processing Techniques 

A final candidate for systematic errors resides in the 
post-processing techniques performed on the spectra. 
The possibilities include incorrect extinction corrections, 
galaxy-contamination corrections, and creation of the 
composite spectra. 

There are two separate extinction corrections applied 
to the spectra. The first is a Mil ky Way extinct i on cor - 
rection, using the dust maps of ISchleeel etall (fl998h . 
The extinction measurements from these maps are typ- 
ically quite good. Furthermore, the Keck/SDSS SNe 
tend to come from low-extinction regions (E(B — V) < 
0.05 mag). It is therefore unlikely that the Milky Way 
extinction corrections are sufficiently erroneous to signif- 
icantly influence the measured flux. 

The second extinction correction is to account for the 
line-of-sight extinction due to host-galaxy dust. The ex- 
tinction has been estimated during the MLCS fitting. For 
our Keck/SDSS sample, the estimates are typically low 
(median value for the Nominal sample is Ay = 0.04 mag 
and the maximum value is Ay — 0.43 mag). How- 
ever, the reddening is determined by comparing the high- 
rcdshift SN light curves to template light curves derived 
from a low-redshift training set. Consequently, if there 
is a difference between the low and high-redshift sam- 
ples, this estimate may be incorrect. On the other hand, 
since the extinction estimates are low, we are making 
small corrections to the spectra, and it is difficult to ex- 
plain how these small corrections could create the large 
UV e xcess in the Keck/S DSS composite spectrum. Fi- 
nally, iFolev et ail (|2008aT ) showed that composite spec- 
tra formed from only low-extinction, uncorrected spec- 
tra and composite spectra formed from both low and 
high-extinction, extinction-corrected spectra are consis- 
tent with each other. 

The light curves of the Keck/SDSS sample were fit with 
Rv = 1.9, the best-fit value found for the full SDSS-II 
samplcF^I. However, this val ue of Ry is low comp ared to 
that seen in the Milky Way (jSchlegel et al.lll998l) . which 
may be partially th e result of mixing SN e la with differ- 
ent intrinsic colors (jFolev fc Kasenll2011h . We have per- 
formed a separate set of light-curve fits with Ry = 3.1. 
Since the value of Ay,host is low for the Keck/SDSS SNe, 
changing the value of Ry did not dramatically alter the 
values of Ay,host- In Figure [TBI we show the Keck/SDSS 
composite spectrum presented in Figure [5] (created from 
the Nominal sample with Ry = 1.9 for all light-curve fits 
and our extinction corrections) and a composite spec- 
trum created from the Nominal sample with Ry =3.1 
for all light-curve fits and extinction corrections. These 

25 This value of Ry was determined by c omparing the SDSS-I I 
photometry to simulations; sec Section 7.2 of Kcsslcr et al. (2009a). 
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Fig. 16. — Composite spectra from Keck/SDSS SNe in the Nomi- 
nal sample with Ry = 1.9 (black curve) and Ry = 3.1 (blue curve) 
for all light-curve fits and extinction corrections. The spectra are 
nearly indistinguishable, indicating that the choice of Ry does not 
affect the UV excess. 

two spectra are created from the same sample of SNe, 
and the only difference is in the value of Ry and de- 
rived Ay,host- Although the different value for Ry does 
change the SED of the composite spectrum slightly, it is 
not a sufficiently large difference to compensate for the 
UV excess. 

As shown in Figure 1191 the galaxy-contamination cor- 
rection appears to be quite good for the observer-frame 
gri bands. These measurements do not include the 
observer-frame near-UV region, where the rest-frame 
UV is observed. However, the galaxy spectra arc de- 
termined from ugriz photometry, which does include 
the observer-frame near-UV region, and other than at 
0.27 < z < 0.33, the reconstruc tion of galaxy SEDs 
is very good (iBlanton et al.|[2003l ). Finally, as seen in 
Figure U the composite spectra created from the full 
Keck/SDSS sample presented in this paper and from the 
sample of SNe with little or no galaxy contamination are 
consistent. It is thus unlikely that incorrect galaxy con- 
tamination causes the UV excess. 

The final possible source of post-processing systematic 
error is related to creating the composite spectra. Since 
the low-redshift and Keck/SDSS composite spectra were 
created in exactly the same way, systematic differences 
should be minimal. The one major possible system- 
atic error associated with creating the composite spectra 
comes from splicing some low-redshift spectra together. 
Unlike with the Keck/SDSS spectra, which all contain a 
large overlapping rest-frame wavelength range, some low- 
redshift spectra are completely disjoint. Accordingly, we 
splice the spectra together using an iterative process, cre- 
ating composite spectra, rescaling spectra, and recreat- 
ing composite spectra until all spectra are added. This 
is especially important for those particular UV spectra 
whose longest wavelengths are ^3200 A. 

For most composite spectra, the overlap between the 
UV and ground-based optical spectra is ~100 A. Conse- 
quently, there may be an offset between the UV and op- 
tical portions of the composite spectrum at ^3200 A (the 
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Fig. 17. — Low-redshift composite spectra with (black curve) and 
without (blue curve) splicing spectra (as described in the text). 
The small differences between the spectra are mainly the result of 
slightly different samples and the iterative scaling of each spectrum 
allowed in the splicing method. The spectra are extremely similar 
for the overlapping wavelengths, indicating that splicing the UV 
spectra to the optical spectra does not affect the flux at the shortest 
optical wavelengths. 



number of spectra for the composite spectra as a function 
of wavelength can be seen in Figure [SJ . We have created 
a composite spectrum without splicing any spectra to- 
gether (using only the SNe that overlap over the same 
1000 A region), and compare in Figure [T7] that spec- 
trum to our composite spectrum consisting of all spec- 
tra. From this figure we see that the method of splicing 
spectra together does not significantly change the com- 
posite spectrum over the wavelength region covered by 
both spectra. Figure [17] also shows that the UV spec- 
tra appear to be spliced properly, creating the correct 
SED. However, we cannot rule out the possibility of er- 
rors at the ends of our spectra combining to incorrectly 
set the scaling of the UV spectra. Regardless of these 
difficulties, there are differences between the Keck/SDSS 
and low-redshift composite spectra at wavelengths longer 
than 3200 A, indicating that there are probably some real 
differences between these samples. With additional UV- 
optical spectra obtained at exactly the same epoch, wc 
may be able to address this further. 

7.2. Supernova Evolution 

No apparent systematic error is creating the UV ex- 
cess for the Keck/SDSS SN la sample. We now investi- 
gate the possible physical reasons for a UV excess in the 
high-redshift SNe la compared to low-redshift SNe la. 
Changes in SNe la with redshift are typically called "evo- 
lution," but it is possible that the sa mple demographics 
are simply shifting wi th redshift (see iFolev et all 2008a 
and ILeibun dgut 200j] for discussions) . We will not dis- 
tinguish between these scenarios here, but we point out 
that if the cause of the difference is related to host-galaxy 
mass (or star-formation rate, etc.), which is compounded 
by the differences between targeted (at low redshift) and 
untargeted (at high redshift) SN searches, this would not 
be considered evolution. 



IK09I found that including and excluding the observed 
u band for nine low-redshift (0.04 < z < 0.09) SDSS- 
II SNe resulted in the same distance modulus offset as 
found for the rest of the SDSS-II sample when including 
and excluding the rest-frame U band. The non-SDSS- 
II low-redshift sample {z w 0.03 for most SNe) did not 
show this offset; if evolution is causing a difference in the 
SEDs, the evolution must occur quickly and recently. 

There are several models for SNe la evolving with red- 
shift as the result of changing progeni t or metallicity (e.g . , 
IHoflich et aL|[l99l ILentz et all l2000l; ISauer et al.|[2008h . 
A major conclusion of these studies is that the optical 
spectra will have relatively small differences, while the 
UV will show a significant difference. Unfortunately, 
there is no consensus as to whether the UV flux will 
increase or decre ase with incre a sing m etallicity. 

The models of IHoflich et al.1 (| 1998( 1 indicate that UV 
excess is ex p ected with higher progenitor metallicity. 
ILentz et"aL1 ([2000), however, showed model spectra 
where UV excess is the result of lower progenitor metal- 
licity. The differences are the result of differin g den sity 
structures (jLentz et al.ll2000HDommguez et al.ll200lD . It 
is beyond the scope of this study to determine which 
model is more accurate. Although other factors such as 
the progenitor C/O ratio may influence the UV flux, the 
difference in UV spectra appears to be broadly consistent 
with models of differing metallicity. 

The slight apparent difference between SN la UV spec- 
tra from SNe with active/low- mass and passive/high- 
mass host galaxies indicates that SNe la with different 
progenitors or delay times may result in different explo- 
sions despite having similar light-curve shapes. It is nec- 
essary to obtain further observations to increase the sam- 
ple size and see if this trend persists. 

There is a clear difference in the relationship between 
TZuv and My for the low-redshift and Keck/SDSS sam- 
ples. It is possible that this is a sign of evolution. How- 
ever, neither the low nor high redshifts have been ad- 
equately sampled in the full My-A plane. We must 
therefore expand these samples to determine if these cor- 
relations persist. However, at this time the relationship 
be tween IZjjv and My found for the low-redshift sample 
by iFFJl clearly should not be used for determining the 
peak luminosity of a high-redshift SN la. 

When comparing the Keck/SDSS and low-redshift 
samples, we have controlled for most observables. In 
particular, the composite spectra shown in Figure[5]have 
the same average phase and light-curve shape. The one 
obvious observable which is vastly different between the 
samples is host-galaxy mass. As we discussed above, the 
average host for an SN in the low-redshift sample is prob- 
ably 10-50 times more massive than that of an SN in 
the Keck/SDSS sample. Given the similarity in other 
parameters, we wonder if such a difference is the domi- 
nant factor in the observed differences. Future samples 
at both low and high redshifts should be able to answer 
this question. 

8. CONCLUSIONS 

Using samples of high-quality low-redshift and 
intermediate-redshift SN la spectra, we have provided 
evidence of differences between the SNe la in these sam- 
ples. By using SN and host-galaxy photometry, we are 
able to properly separate the SN and host-galaxy spectra. 
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Constructing a composite spectrum from 17 of our host- 
galaxy contamination-corrected Keck/SDSS SNe, we see 
that the intermediate-redshift SNe have > 20% more 
UV flux (A < 3400 A)of their low-redshift counterparts. 
This UV excess is also apparent in synthesized broad- 
band colors at high significance. After an exhaustive 
search for systematic errors, we determined that poor 
spectral reductions, flawed post-processing techniques, 
or Malmquist-biased sample selection are unlikely to be 
creating this effect. It is possible, however, that a combi- 
nation of several small systematic effects creates the UV 
excess. Despite having otherwise very similar properties, 
the UV flux of the Keck/SDSS sample does not appear 
to be representative of the larger SDSS-II sample. 

We have also shown that there may be a slight differ- 
ence in the UV spectra of SNe from low-mass and high- 
mass host galaxies (or, alternatively, those that are ac- 
tively forming stars versus those that are not). This dif- 
ference is currently of low statistical significance and re- 
quires larger samples to further improve confidence in the 
result. In particular, SNe with more massive hosts must 
be probed at high redshift. If SNe la originating from dif- 
ferent stellar populations have similar optical light-curve 
shapes but different UV spectra, observing the UV spec- 
tra of SNe la may offer a new window to understanding 
the progenitors of SNe la and could potentially improve 
distance estimates. Several studies have now found a 
significant correlation between host-galaxy mass, mor- 
phology, and star-formation rate with Hubble residuals. 
Detailed analysis of spectra may reveal the physical cause 
of these correlations. 

The UV ratio, IZjjv, has a substantially different corre- 
lation with My in the low-redshift and Keck/SDSS sam- 
ples. Since these samples are relatively small and that 
of Keck/SDSS has little diversity in SNe la, it is still 
unclear whether this difference is inherent to low and 
high-redshift samples or is the result of insufficient sam- 
ples. More observations at both low and high redshifts 
are necessary to make this determination. 

The differences in the UV spectra of low and high- 
redshift SNe la may be the result of SN la evolution. 
Future observations and careful modeling are necessary 
to determine the magnitude and implications of this pos- 
sibility. However, a more reasonable explanation is that 
the low-redshift and high-redshift samples differ in some 
significant way. We have accounted for phase and light- 
curve shape when comparing these samples, but given 
the heterogeneous nature of (and incomplete information 
for) the low-redshift sample, we were unable to account 
for additional potential differences such as host-galaxy 
properties. Since the low-redshift SNe la were discovered 
through surveys which target the most luminous galaxies 
in the nearby Universe and the SDSS-II survey was un- 
targeted, it is likely that the host-galaxy properties for 
these samples are quite different. 

The UV excess seen in the Keck/SDSS sample is not 
seen in the full SDSS-II sample, as inferred from their 
light curve (K09). In fact, the full SDSS-II sample is, on 
average, UV deficient relative to the low-redshift sam- 
ple. There is no clear reason for this difference, and it 
may simply be the result of a relatively small sample (al- 
though the Keck/SDSS sample represents 12% of the full 
|KQ| SDSS-II sample). 

A mismatch in the UV SEDs of SNe la like that seen 



between the Keck/SDSS and the low-redshift samples 
has important implications for the estimated distances 
to all SNe la. Although the Keck/SDSS sample is rel- 
atively large, it still is not fully representative of the 
larger SDSS-II sample. This indicates that there is a 
significant scatter in the UV properties of SNe la (at 
least at intermediate redshift). A simplistic treatment of 
this flux difference suggest that the incorrect normaliza- 
tion of the peak luminosity (an offset in My iPca k for a 
given light-curve shape) and an incorrect determination 
of the host-galaxy extinction (changing the measurement 
of Ay) should both contribute to an underestimate of the 
correct distance for an SN la. Fitting multi-band data 
that cover the rest-frame UV and optical will tend to 
cause bluer SNe la to have a higher measured peak lumi- 
nosity than reality, causing an opposite (and dominant) 
effect. This effect is constrained to the (near) UV of 
an SN la SED, and by ignoring this region, one should 
obtain an unbiased (by this effect) distance estimate. 

A light-curve fitter trained on a low-redshift sample 
(MLCS2k2) and applied to high-redshift samples would 
produce significantly different distance measurements if 
the rest-frame U band were included or excluded. This 
was not seen with a light-curve fitter which was trained 
on the full (low and high-redshift) sample (SALT2). If 
the cause of the [/-band anomaly was a mismatch in the 
UV SEDs for the low and high-redshift samples, SALT2 
would have a model SED close to both the low and high- 
redshift samples and differences would be minimized. 

For the James Webb Space Telescope (JWST), Eu- 
clid, and the future potential Wide-Field InfraRed Sur- 
vey Telescope (WFIRST) mission, which will observe 
SNe in the near-infrared, there should be a relatively 
small systematic color bias to a redshift of ~2.5. Since 
current designs do not plan to observe SNe la beyond 
that redshift, differences seen in our sample should have 
a relatively small effect on JWST, Euclid, or WFIRST 
measurements. 
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APPENDIX 
GALAXY SUBTRACTION 
Photometry Matching and Spectral Warping 

The first method to correct for galaxy contamin ation makes no as sumptions about the quality of our spectropho- 
tometry and is similar to the method presented by lEllis et al.l (|2008( ). We perform extractions of our SNe la with a 
fixed aperture and background regions which avoid any host-galaxy light. Our resulting spectra are a combination of 
SN and host galaxy with no attempt to remove host-galaxy contamination. 

As a byproduct of the SMP photometry, we have th e broad-band p hotometry for the host galaxy at the position of 
the SN. Using multicolor light-curve shape (MLCS2k2; Jha et al. 2007) template light curves, we are able to interpolate 
the SMP photometry (independently in each band) to determine the broad-band photometry for the SN at the time 
the spectrum was obtained. Using the photometry of both the SN and the background at the position of the SN, we 
warp our spectrum to match the measured photometry of the SN and background. We do the warping using both a 
spline and linear fit with three pivot points at the flux-weighted centers of the gri filters. Since the SEDs are slightly 
different for each object, we allow the pivot points to shift in wavelength to match the true center of the filter for each 
spectrum. We can then subtract the reconstructed galaxy SED (properly scaled to the galaxy photometry) from the 
spectrum, yielding a spectrum which matches the photometry of the SN. 

Since a spline will exactly fit the measured photometry, we perform a Monte Carlo simulation to determine if the 
spline fitting is producing reasonable residual spectra. To do this, we randomly choose photometry for each object 
assuming Gaussian distributions for the photometry with a mean corresponding to the measured photometry and a 
width corresponding to the uncertainty in the measurements. Since the uncertainties in the gri bands are of similar 
magnitude, the average of the residuals from the Monte Carlo simulations is typically very similar to the spline using 
only the measured photometry. 

We find that the photometry is always well fit by a linear function of wavelength resulting in % 2 /dof >C 1, indicating 
that for many SNe this method is ovcrfitting the data (i.e., no galaxy subtraction is required for some SNe). For most 
cases the slope of the line is < la (and in all cases the slope is < 3tr) from 0, again indicating little to no correction 
required. For our entire sample, there is no trend in the slope for a given night, standard star, or slit position angle, 
with nearly equal numbers of SNe with positive and negative slopes, suggesting that there is no systematic problem 
affecting our galaxy subtraction. 

Color Matching 

The second method initially assumes that our spectrophotometry is well calibrated (i.e., errors from the spectropho- 
tometry do not dominate the error in comparing the spectrum to the photometric colors). The flux calibration is 
verified both by the photometry matching above and also after correcting for host-galaxy contamination. For this 
method, we attempt to remove as much galaxy contamination as possible when extracting our spectra. Although 
removing galaxy light at the point of extraction prevents matching to the absolute scaling of photometry, it tries to 
reduce the introduction of errors by subtracting a potentially imprecise galaxy SED from a spectrum. Since the galaxy 
photometry determines the reconstructed galaxy SED, if our spectrophotometry is well calibrated, then subtracting 
the reconstructed galaxy SED from the observed spectrum such that the colors of the residual spectrum matches those 
of the SN will result in an SN-only spectrum. 

We can demonstrate this mathematically as follows. In general, an observed SN spectrum is defined by 

/s P cc = A(/ SN + B/ ga i), (Al) 

where / sp cc, /sn, and / ga i are the vectors of fluxes in the observed spectrum, SN spectrum, and galaxy spectrum, 
respectively, and A and B are normalization factors. One can think of A as normalizing the spectrum in an absolute 
sense to account for slit losses, clouds, and other effects. The parameter B can be thought of as a color normalization, 
where it is adjusted until / spe c and /sn + Bf ga \ have the same color. 

As a byproduct of the SMP photometry, we have p ga i, the broad-band photometry (in flux units) for the host galaxy 
at the position of the SN. Using multicolor light-curve shape (MLCS2k2; |jha et al.ll2007r i template light curves, we are 
able to interpolate the SMP photometry (independently in each band) to determine psn , the broad-band photometry 
(in flux units) for the SN at the time of the spectrum. 
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We can define the function which translates spectra to synthetic broad-band photometry as P, where P(/sn) = Psn 
and P(f g ed) = P g ed- Note that we impose the first relationship, while the second relationship is required by our method 
of determining / ga i . 

From our spectrum, we are able to determine P(/ sp oc), the broad-band synthetic photometry (in flux units) of the 
spectrum, which includes both SN and galaxy light. These vectors obey the equation 

P(/ SP ec) = A few + B Pgal ) . (A2) 

For this equation to be valid, we make two assumptions. The first, already noted above, is that our spectra are well 
flux calibrated. The second assumption is that B, the relative fraction of the galaxy and SN light in the observed 
spectrum, does not vary strongly with wavelength. We will return to both of these assumptions below. 

Having gri photometry for the galaxy and the SN, as well as the observed spectrum and the reconstructed galaxy 
spectrum, we are able to determine the factors A and B. We can then use these factors to appropriately remove the 
reconstructed galaxy spectrum from the SN spectrum, 

/SN = A / S p ec - P/gal- (A3) 

We have applied this technique to all SN spectra. Figure Q] displays two of our observed spectra, the reconstructed 
host-galaxy spectra, and the galaxy-contamination corrected SN spectra (see the online version for figures showing all 
spectra). 

Galaxy Subtraction Performance 

We have applied two methods to properly correct for host-galaxy contamination and to recover the SN SEDs for 
our sample of SNe la. One method attempts to match spectra to the measured brightnesses, while the other strives 
to match to the measured colors. Neither technique is superior for all cases, but we will discuss their performance for 
our data below. 

The greatest advantage of the photometry-matching method is that it can minimize errors from poor spectropho- 
tometry. By definition, the spline-fitting version yields SN SEDs that exactly match the SN photometry, while the 
line-fitting version will never increase the residuals. However, this method has two main drawbacks. First, it makes 
no attempt to reduce galaxy contamination during extraction. Second, the warping of the spectrum is somewhat 
arbitrary and must be extrapolated to wavelengths beyond those covered by our broad-band photometry. 

In a situation where a galaxy has a constant SED at all positions and a smoothly varying surface-brightness profile, 
it is simple to remove the galaxy contamination in a spectrum during extraction. By measuring the galaxy light on 
either side of the position of the SN in the two-dimensional spectrum and interpolating over the position of the SN, one 
can cleanly remove all galaxy contamination. If, however, one makes no attempt to remove the galaxy contamination 
during extraction, but rather relies on a reconstruction of the galaxy SED from broad-band photometry, one will 
almost certainly introduce additional errors in the spectrum. 

If our spectrophotometry is incorrect for a given object, the correct function to warp the spectrum to match the 
photometry can never be perfectly determined. For our sample, we have assumed that this function is either a spline 
or a line with anchors at the center of the gri bandpasses with fluxes equal to those measured from the gri photometry. 
Although the corrections were typically small, deviations from the fit are possible for wavelengths between the anchor 
points. Additionally, any wavelengths beyond the anchor points must be extrapolated. Since our wavelength coverage 
(^3200-9200 A) is much larger than the span from the bluest to the reddest anchor points (approximately 4800 and 
7500 A for g and i, respectively), we must either extrapolate for a significant portion of our spectra or trim our 
spectra. For the spline fitting, the implied corrections in the extrapolation region can be large, while the corrections 
from the linear fitting are more realistic. For the spline-fitting version, we have chosen to extrapolate to the blue edge 
of the g band and the red edge of the i band, 3950 A and 8220 A, respectively. For these wavelengths, we have some 
information contributing to the photometry, while also limiting our extrapolation. For the linear-fitting version, we 
extrapolate further, to 3200 A and 9000 A. 

The one key advantage of the color-matching method is that galaxy subtraction is performed during the extraction 
of the spectra. This will reduce any error associated with having an incorrect galaxy SED. Again, this method requires 
two assumptions: the spectrophotometry is not the dominant source of error and the PSF of our object does not vary 
wildly with wavelength in a way where our parameter B becomes strongly dependent on wavelength. 

We can examine our spectrophotometry by comparing synthetic broad-band colors of our spectra to colors from 
the photometry. In particular, we are most interested in the relative spectrophotometry of our spectra without any 
attempt to remove galaxy contamination during spectral extraction. These spectra can be directly compared to the 
photometry of the SN and galaxy at the position of the SN. 

In Figure 1181 we present the g — r and r — i colors one would find by measuring the fluxes in an aperture centered 
on the SNe at the time of our spectra with no host-galaxy subtraction in our images, as well as the synthetic colors 
of our spectra without any attempt to remove galaxy contamination during spectral extraction. These measurements 
are consistent with x 2 /dof = 21.15/21, showing that our spectrophotometry is very good over the bands where we can 
measure it for these SNe. 

A further test of our spectrophotometry is to examine the synthetic photometry of our SN spectra after the galaxy 
contamination has been removed. In the case of the photometry-matching technique, the residuals are by design 
small. For the color-matching method, on the other hand, the residuals are not guaranteed to be small. If our 
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Fig. 18. — (Left) Synthetic g — r colors of our SN spectra with no attempt to remove galaxy contamination compared to the photometric 
g — r colors determined by adding the photometry of the SN at the time of the spectrum to that of the galaxy at the position of the SN. 
Overplotted is the line indicating equivalent photometric and synthetic colors. (Right) Synthetic r — i colors of our SN spectra with no 
attempt to remove galaxy contamination compared to the photometric r — i colors determined by adding the photometry of the SN at the 
time of the spectrum to that of the galaxy at the position of the SN. Overplotted is the line indicating equivalent photometric and synthetic 
colors. The assumption that our spectrophotometry is correct yields x 2 /dof = 21.15/21. 



spectrophotometry (or alternatively, our galaxy SEDs) were incorrect, we could obtain a minimal residual in the colors 
that are still far from the measured values for the SN. In Figure [T9l we present the ratio of the synthetic fluxes recovered 
after galaxy subtraction to the measured photometric fluxes for each SN in our sample. Ignoring SNe la with low S /N 
or failed galaxy-template determination, the average of the residuals is 1.03, 0.97, 1.01, with a scatter in the residuals 
of 5.4%, 4.5%, and 5.0% in the g, r, and i bands, respectively. If we also remove the SNe with 0.27 < z < 0.33, 
the average of the residuals is 1.03, 0.97, 1.01, with a scatter in the residuals of 5.1%, 3.7%, and 4.2% in the g, r, 
and i bands, respectively. These un certainties are consist ent with the errors expected from the spectroph otometry, 
photometry, and galaxy subtraction (Matheson ct al. 2008; Silverman et al., submitted; Blanton et al.ll2003l ). For the 
photometry-matched, linearly corrected spectra including (excluding) SNe with 0.27 < z < 0.33, the average of the 
residuals is 1.06, 1.06, and 1.08 (1.04, 1.03, and 1.07) with scatter in the residuals of 6.0%, 9.0%, and 11.8% (3.0%, 
6.6%, and 9.4%) for the g, r, and i bands, respectively. 

From Figure [191 we can also infer that our parameter B does not vary significantly with wavelength. If, for instance, 
B were a monotonically increasing function of wavelength, but we were not accounting for it, we would expect to 
oversubtract our galaxy SED in the g band and undersubtract our galaxy SED in the r band. Since the average values 
of the ratio of the synthetic to photometric flux are consistent with 1.0 in all bands, we conclude that B is not strongly 
dependent on wavelength. 

From our tests, we find that both the photometry-matching and color-matching methods produce excellent galaxy- 
contamination corrected SN spectra over the gri bands. Because of the extrapolation of the spline function, we do 
not believe our spline-fit photometry-matched spectra can be extended blueward of the g band or redward of the 
i band; similarly, the linear-fit photometry-matched spectra are limited to 3200-9000 A. Considering the excellent 
spectrophotometry of our spectra for the gri bands and our experience with t/-band flux calibration for many low- 
redshift SNe la (Silverman et al., submitted), the color- matching technique should yield properly flux-calibrated spectra 
redward of ^3200 A. We have performed our analysis using both the photometry-matched and color-matched spectra, 
yielding similar results. We will present some results from both methods, but will focus on the color-matched spectra. 

Galaxy Subtraction Failures 

An SN fails the galaxy subtraction when the best-fit galaxy SED has larger flux values than the galaxy-contaminated 
SN spectrum for a significant portion of the spectrum. This can occur if there is incorrect galaxy or SN photometry or 
if the SN spectrum is not properly flux calibrated. Additionally, for SNe with 0.27 < z < 0.33, where incorrect galaxy 
SEDs are possible, the reconstructed SED may be unphysical. 

There are four SNe la for which the galaxy subtraction fails. Three of these (SNe 16644, 16789, and 20829) have 
redshifts of ~0.3 and significant galaxy contamination. The unphysical SEDs of reconstructed galaxy spectra at 
z » 0.3 are the obvious reason for a failure of the fit. Despite this, by constraining a heavily smoothed version of the 
residual SN spectrum to have non-negative fluxes at all places besides the ends and at emission lines (which may be 
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Fig. 19. — Ratio of synthetically created fluxes to measured photometric fluxes for the Nominal sample. The black crosses, blue diamonds, 
and red triangles correspond to SNe with no, low (< 8%), and high (> 8%) host-galaxy contamination, respectively. SNe where the galaxy 
subtraction did not fail, but are at redshifts 0.27 < z < 0.33, have been marked with an additional red square. The dotted lines represent 
an estimate of the expected error from a combination of photometry, spectral-flux calibration, and galaxy-subtraction errors. 

oversubtracted), we can achieve a rough, but probably oversubtracted, SED of the SN. Although this is not a rigorous 
treatment of the galaxy contamination, it could be useful for visualization. 

SN 7512 also fails the galaxy-subtraction procedure. It has relatively low S/N spectra, making galaxy subtraction 
less precise. The four SNe with failed galaxy subtraction have been removed from our sample for any analysis, leaving 
a sample of 17 SNe la that we use for our analysis. Additional spectroscopy of the SN hosts could improve the SEDs of 
the SNe with 0.27 < z < 0.33, but the faintness of the unobserved host galaxies makes this task difficult with current 
facilities. 

COMPARISON TO SNLS Z a 0.5 COMPOSITE SPECTRUM 

Two recent studies generated composite spectra of high-redshift SNe la . iFolev et al.1 (l2008afl presented composite 
spectra generated from high-redshift SNe la from the ESSENCE survey (|Miknaitis et al.l 120071 ). Unfortunately, the 
lack of multi-color information prevented galaxy-contamination correction similar to that impleme nted in this study . 
As a result, there was little information on changes to the SN SED. However, a second study by I Ellis et al.l (|2008l ) 
of SNLS SNe la did present composite spectra after correcting for the galaxy contamination. The two spectra (one 
"early" spectrum with t < —4 days and one "maximum" spectrum with —4 < t < 4 days; both with (z) ~ 0.5) can 
provide an interesting comparison to our composite spectra. 

In Figure [20l we illustrate the Keck/SDSS, SNLS, and low-redshift maximum-light composite spectra. For this 
comparison, we scaled the spectra to the SNLS spectrum over its entire wavelength range. The SNLS spectrum is 
consistent with both the Keck/SDSS and low-redshift spectra over the wavelength range sampled. The early and 
maximum-light SNLS spectra both have a median stretch of 1.06, which corresponds to A = —0.23. Therefore, the 
SNLS SNe have broader light curves (and are correspondingly more luminous SNe la) than those of the Keck/SDSS or 
low-redshift samples. The UV SEDs of SNe la are strongly dependent on their light-curve shape (see Section l5~3|) . and 
since we do not know the average light-curve shape of the composite spectrum (the individual spectra have different 
weights and therefore contribute different amounts to the composite spectrum) , the comparison is of less utility than 
that of the low-redshift and Keck/SDSS composite spectra. Despite these differences, it is worth noting that the SNLS 
spectrum has less flux than the Keck/SDSS spectrum (but still within the error region) and slightly more flux than 
the low-redshift spectrum at most wavelengths (again, still within the error region). Larger differences may exist in 
these samples, but since the SNLS spectrum does not probe much of the optical region (stopping blueward of 5400 A), 
the lever arm of scaling the spectra to match at long wavelengths is unavailable for this analysis. As seen with the 
low-redshift and Keck/SDSS spectra in FigureHOl matching the spectra at only these blue wavelengths can cause large 
divergence at redder wavelengths. 

In Figure [2X1 we present a comparison of the SNLS "Early" composite spectrum with a low-redshift (t) = —6.6 days 
composite spectrum from IFolev et al.l ([2008a). Over the restricted wavelength range, the two spectra appear to be 
consistent. There is a large deviation at ^3100 A, with the SNLS spectrum having excess flux compared to the 
low-redshift spectrum, but the uncertainties in both the SNLS and low-redshift spectra are large enough for the 
spectra to be consistent. Similar to the maximum- light spectra, these samples likely have different light-curve shape 
characteristics, and additional information is necessary to make definitive conclusions. 



28 



Foley et al. 



Oi 0.4 



SDSS 

Low-Redshift 
SNLS 




3000 4000 5000 6000 

Rest Wavelength (A) 



7000 



Fig. 20. — (Top panel): Co mposite spectrum created from our Keck/SDSS sample (black curve) compared to the SNLS maximum- 
light composite spectrum from [Ellis et al.l (:2008) (red curve), and the maximum-light low-redshift composite spectrum (blue curve). The 
grey, yellow, and light-blue regions are the la boot-strap sampling errors for the Keck/SDSS, SNLS, and low-redshift composite spectra, 
respectively. (Bottom panel): The yellow region is the la boot-strap sampling region for the SNLS composite spectrum. The black and 
blue curves are the residuals of the SNLS and the Keck/SDSS (Keck/SDSS minus SNLS) and low-redshift (low redshift minus SNLS) 
composite spectra, respectively. The grey and light-blue regions are the residuals of the SNLS composite spectrum and the Keck/SDSS 
and low-redshift la boot-strap sampling regions, respectively. 




£ -°- 2 



3000 



3500 



4000 4500 
Rest Wavelength (A) 



5000 



Fig. 21. — (Top panel): SNLS "Early" composite sp ectrum created from [Ellis et al. (2008) (red curve; t < —4 days) compared to the pre- 
maximum-light low-redshift composite spectrum from Foley et al. (2008a) (blue curve, (t) = —6.6 days). The yellow and light-blue regions 
are the lcr boot-strap sampling errors for the SNLS and low-redshift composite spectra, respectively. (Bottom panel): The yellow region 
is the lcr boot-strap sampling region for the SNLS composite spectrum. The blue curve is the residual of the SNLS and the low-redshift 
composite spectra. The light-blue region is the residual of the SNLS composite spectrum and the low-redshift la boot-strap sampling 
regions. 



The expectation is that any physical change in SN la spectra with redshift would be monotonic with redshift. With 
that in mind, the fact that the SNLS composite spectrum is slightly above the low-redshift composite spectrum and 
slightly below the Keck/SDSS composite spectrum is puzzling. Perhaps the Keck/SDSS sample is biased in a way 
we have not been able to detect, or our assumption of monotonic spectral behavior is incorrect. However, it is worth 
noting that when splitting the Keck/SDSS sample by redshift, the higher-redshift composite spectrum has less UV flux 
than the lower-redshift composite spectrum. Although this difference is not significant given our boot-strap sampling 
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errors, it shows the same trend as with the SNLS spectra. 

Without combining the individual SNLS SN la spectra in the same manner as t he Keck/SDSS SN la spectra, it is 
difficult to quantify the similarities and differences. In particular. lEllis et al.l (|2008[ ) warped the spectra to match their 
photometry before correcting for host-galaxy contamination. The reasoning for this is that differential slit losses and 
atmospheric extinction cause errors in the spectrophotometry. We have show n that our reduction m ethods reproduce 
B — V colors with a low dispersion (cr = 0.063 mag) for low-redshift SNe la (jMatheson et al.l 12008t h and we see that 
there is a low scatter in the synthesized photometry (Figure [T9l of our sample of objects with no discernible host galaxy 
in the images (despite there possibly being some underlying host contributing slightly to the spectrum). Hence, we 
conclude that our spectrophotometry is reasonable an d has uncer t aintie s similar to those associated with host-galaxy 
contamination. Furthermore, the spectral warping of lEllis et al.l (|2008l ) extrapolates into the rest-frame UV for the 
lower-redshift (z < 0.3) objects. 



